AD— All  1  422  ILLINOIS  1>NIV  AT  URBANA  ENGINEERING-PSYCHOLOGY  RESEAR—ETC  F/G  5/10 

THE  PRODUCTION  OF  SPECIFIED  ELECTROCORTICAL  ACTIVITY  AS  A  MEASU — ETC(U) 
FEB  B1  0  B  BER INGER »  M  G  COLES  DACA88-60-C-0003 

UNCLASSIFIED  EPL-01-2/USACERL-81-1  NL 


f 


/  \ 


Report  No.  EPL-81-2/USACERL-81-1 


The  Production  of  Specified 
Electrocortical  Activity 
as  a 

Measurable  Task 


TECHNICAL  REPORT 


Prepared  for:  The  U.  S.  Army  Construction  Engineering 
Research  Laboratory  and  the  Air  Force 
Office  of  Scientific  Research 


Prepared  by: 


Dennis  B.  Beringer 

Engineering  Psychology  Research  Laboratory 
University  of  Illinois  at  Urbana-Champaign 

Michael  6.  H.  Coles 

Department  of  Psychology 

University  of  Illinois  at  Urbana-Champaign 


opprov^ 

for  pubWe  r»l  <wd  saJe:  U» 

distribution  l»  wibadtodt 


4 


03 


01  019 


82 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entmrod) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

EPL-81 -2/USACERL-81  -1  /)  h-  f\  XL 

3.  RECIPIENT'S  CATALOG  NUMBER 
/  / 

/■ -  W 

4.  TITLE  (mid  Subttllo) 

THE  PRODUCTION  OF  SPECIFIED  ELECTROCORTICAL 

ACTIVITY  AS  A  MEASURABLE  TASK 

5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

TECHNICAL 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORS) 

Dennis  B.  Beringer 

Michael  G.  H.  Coles 

B.  CONTRACT  OR  GFfcANT  NUMBERfm) 

WA ; 

DACA-88 -80 - C -000 3 

!•,  PERFORMING  ORGANIZATION  .name  anix  address.  _  _ 

The  University  of  Illinois,  Dept,  of  Indus.  Engr 
Engineering  Psychology  Laboratory 

Champaign,  Illinois  61820 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

]  ]  /  A  -LiO-i/i  ATI  4, 

(/:  research  lab 

Interstate  Research  Park,  Champaign,  IL  61820 

Air  Force  Office  of  Scientific  Research  (NL) 

Boiling  Air  Force  Base.  DC.  20332 

12.  REPORT  DATE 

February  1981 

12.  NUMBER  OF  PAGES 

60  fK 

14.  MONITORING  AGENCY  NAME  «  ADDRESSf//  tUUmtont  fnm  Controlling  out  do) 

Same 

15.  SECURITY  CLASS,  (of  Mm  roport) 

UNCLASSIFIED 

iBo,  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

1«.  DISTRIBUTION  STATEMENT  (of  ihl •  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  tho  mbmtrmct  mntormd  In  Block  20,  ft  difforont  fnm  Roport) 


IS.  supplementary  NOTES 


I  IS.  KEY  WORDS  (Continue  on  eoworoo  mldo  If  nocooeory  and  identify  by  block  number ) 

Electrocortical  Activity  Feedback  Utilization 

Productivity  Biofeedback 

Alpha  Production  Micropressor  Systems 

Cognitive  Models 

M.  AWTWACT mwiw  lO  Ifiitw— y  wilrfrtiiy  ty  Mort  wiirj  "  ~  * « 

^  Global  electrocortical  activity  (alpha)  was  examined  as  a  possible  Index 
of  subtle  environmental  effects  upon  productivity  in  cognitive  tasks.  A  new 
model  of  biofeedback  was  proposed  to  account  for  the  various  effects  described 
In  the  literature.  Although  consistent  with  previous  findings,  the  model 
received  no  addltibnal  confirmation  from  the  current  study.  It  was  concluded 
that  the  measure  was  overly  sensitive  and  unstable.  Other  multivariate 
measurement, schemes  were  suggested.  ' 


_ jJJJJlr 

SECURITY  CLASSIFICATION  OP  THIS  PAGE  fSE»  Dmto  Stored) 


DO  |  jam  »  1473L  CDfTIOM  O T  1  MOV  SI  tft  OBSOLETE 


I 


FOREWORD 


This  study  was  jointly  funded  by  the  Office,  Chief  of  Engineers 
under  project  4A161102AT23,  and  the  Air  Force  Office  of  Scientific  Re¬ 
search  under  Program  Element  61102F,  Project/Task  2313/A3.  The  study 
was  performed  under  contract  by  the  Engineering  Psychology  Laboratory 
of  the  Department  of  Industrial  Engineering,  University  of  Illinois  at 
Urbana-Champalgn.  The  contract  was  monitored  by  the  Facility  Systems 
Division,  U.  S.  Army  Construction  Engineering  Research  Laboratory, 
Champaign,  IL. 


i 


l 


ii 


TABLE  OF  CONTENTS 


INTRODUCTION 


Models  of  Biofeedback . 2 

Mulholland's  systems  approach  .  2 

Schwartz's  systems  approach  .  4 

A  new  model  of  biofeedback  .  4 

Mediation  .  6 

Hypotheses  .  7 


PHASE 


Hardware  Devices  .  9 

Software  systems  .  11 

Method . 11 

Results  and  Discussion  .  13 

PHASE  2 . 18 

Method . 19 

Results  and  Discussion  .  20 

CONCLUSIONS . 22 

REFERENCES . 23 

APPENDIX  A:  SCHEMATIC  DIAGRAMS  .  24 

APPENDIX  B:  SOFTWARE  LISTINGS  .  30 

APPENDIX  C:  ANOVA  SUMMARY  TABLES  .  47 

APPENDIX  D:  CORRELATIONS  AND  REGRESSION  .  57 


i 


LIST  OF  TABLES 


Table  Page 


1  Summary  of  analysis  of  variance  of  %  alpha  categorical  by  training 

exposure . 48 

2  Summary  of  analysis  of  variance  of  log  cycle  frequency  categorical 

by  training  exposure  .  49 

3  Summary  of  analysis  of  variance  of  %  alpha  categorical  by  feedback 

condition  and  eye  condition . 50 

4  Summary  of  analysis  of  variance  of  baseline-scaled  %  alpha  categorical 

by  feedback  condition  and  eye  condition  . ' .  51 

5  Summary  of  analysis  of  variance  of  log  baseline-scaled  %  alpha  by 

feedback  condition  and  eye  condition  .  52 

6  Summary  of  analysis  of  variance  of  reciprocal  %  alpha  by  feedback 

mode  and  feedback  polarity . 53 

7  Summary  of  analysis  of  variance  of  arcsine  cycle  frequency  by  feedback 

mode  and  feedback  polarity . 54 

8  Summary  of  analysis  of  variance  of  square-root  baseline-scaled  %  alpha 

by  feedback  mode  and  feedback  polarity . 55 

9  Summary  of  analysis  of  variance  of  arcsine  baseline-scaled  cycle  fre¬ 
quency  by  feedback  mode  and  feedback  polarity . 56 


iv 


LIST  OF  FIGURES 

Figure  Page 

1  Feedback  system  with  feedforward  and  feedback  (Mulholland,  1977)  .  .  3 

2  Schwartz's  model  for  biofeedback  (Schwartz,  1979)  .  3 

3  The  proposed  model  of  biofeedback  .  .  5 

4  Data  collection  system  hardware  configuration  .  10 

5  Time  in  alpha  by  time  on  task:  Superposed  trials  (A),  mean  and  s.d. 

(B),  and  hypothetical  data  (C) . 14 

6  Retransformed  condition  means:  Baseline-scaled  %  alpha  by  feedback 

and  eye  conditions . 17 

7  Schematic  diagram:  Analog-to-analog  interface,  channel  1  25 

8  Schematic  diagram:  Analog-to-analog  interface,  channel  2  26 

9  Analog-to-analog  interface  input/output  characteristics  .  27 

10  Schematic  diagram:  Digital  stimulus  controller  .  28 

11  Schematic  diagram:  Digital  tone  generator  .  29 

12  Correlation  and  regression:  scatter  plot  of  %  alpha  by  cycle  fre¬ 
quency  across  all  trials  and  subjects . 58 

13  Correlation  and  regression:  scatter  plot  of  %  alpha  by  cycle  fre¬ 
quency  across  all  subjects  for  eyes-closed  trials  .  59 


« 


14  Correlation  and  regression:  scatter  plot  of  %  alpha  by  cycle  fre¬ 
quency  across  all  subjects  for  eyes-open  trials  . 


60 


INTRODUCTION 


During  the  last  decade,  considerable  attention  has  been  paid  to  the  possible 
application  of  EEG  Biofetdback  Training  in  a  number  of  different  situations. 
However,  research  to  date  has  indicated  that  biofeedback  has  few,  if  any, 
beneficial  effects.  There  have  been  two  major  areas  in  which  the  search  for 
possible  application  has  been  pursued.  In  the  clinical  area,  it  was  believed 
that  the  production  of  alpha  activity  through  biofeedback  techniques  would  pro¬ 
duce  a  state  antithetical  to  many  clinical  disorders  such  as  anxiety  and 
depression.  However,  as  Sterman  (1977)  indicates,  the  promise  of  biofeedback 
in  this  area  has  not  been  realized.  A  second  area  of  disappointment  has  been 
that  of  performance  enhancement.  Lawrence  and  Johnson  (1977)  reviewed  the 
research  in  this  area  and  concluded  that  "enhanced  alpha  activity  does  not 
prevent  sleep  loss  effects  or  substitute  for  sleep  ...  is  not  related  to 
memory  or  choice  reaction  time  performance  .  .  .  does  not  provide  a  recupera¬ 
tive  break  period  .  .  .  and  is  incompatible  with  cognitive  tasks  requiring  any 
degree  of  effort"  (pg.  166).  Thus,  the  promise  of  biofeedback  as  a  panacea  in 
clinical  and  human  performance  applications  has  not  been  realized. 

The  present  studies  take  a  different  approach  to  the  potential  application 
of  biofeedback.  Instead  of  seeking  beneficial  effects  from  feedback  procedures, 
the  approach  seeks  to  evaluate  whether  performance  on  a  biofeedback  task  can  be 
used  as  a  measure  of  the  effects  of  work  environment.  Briefly  stated,  the 
research  is  predicated  on  the  hypothesis  that  the  biofeedback  task  involves  the 
allocation  of  resources  or  capacity.  External  events  (environmental  stimuli) 
or  internal  events  (cognitive  processes)  can  interfere  with  the  biofeedback 
task  since  they  draw  on  resources  involved  in  the  task.  Thus,  performance  on 
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the  biofeedback  task  can  be  used  as  an  index  of  the  degree  to  which  environmental 
events  and  cognitive  events  call  upon  the  same  processes  that  are  involved  in 
the  task.  Before  expanding  on  the  theoretical  basis  for  these  predictions,  we 
need  to  consider  various  models  of  biofeedback. 

tfadgla  sL  mafegtoagk 

Mulholland's  systems  approach .  The  feedback  system  described  by  Mulholland 
(1977)  is  based  on  a  generalized  feedback  control  system  (Figure  1).  Empirically, 
it  has  been  tested  using  biofeedback  of  EEG  with  a  feedback  stimulus  being 
one  which  itself  influences  the  EEG.  In  most  studies,  a  visual  feedback 
stimulus  is  used  and  occipital  EEG  recorded.  The  visual  stimulus  "naturally" 
results  in  an  attenuation  of  alpha  activity  (e.g.,  Gale,  Dunkin,  and  Coles,  1969). 
Now,  under  eyes-open  conditions,  resting  EEG  shows  a  cycling  between  alpha  and 
no-alpha  states.  Thus,  if  the  onset  of  the  visual  stimulus  is  linked  to  the 
production  of  alpha,  the  natural  cycling  between  alpha  and  no-alpha  states 
will  be  influenced.  Mulholland  has  shown  that  variability  in  cycling  is  re¬ 
duced  by  the  provision  of  this  type  of  feedback  system  and,  from  a  series  of 
studies,  has  concluded  that  the  systems  model  shown  in  Figure  1  is  a  reason¬ 
able  way  to  approach  biofeedback  phenomena. 

One  important  point  should  be  made  about  Mulholland's  model.  The  subject 
is  not  required  to  (a)  be  aware  of  the  contingencies  of  the  feedback  or  (b) 
consciously  try  to  influence  the  feedback.  In  this  respect,  Mulholland's 
research  is  not  directly  relevant  to  the  present  studies  where  the  subjects 
are  required  to  invent  some  cognitive  (or  other)  resources  in  the  generation 
of  alpha.  However,  the  research  does  point  to  the  importance  of  the  feedback 
stimulus  as  a  source,  not  only  of  Information  for  the  subject,  but  also  as  an 
event  which  exerts  a  "natural",  unconditioned  effect  on  the  EEG. 


DISTURBANCE 


Figure  1.  Feedback  system  with  feedforward  and  feedback  (Mulholland,  1977). 


Figure  2.  Schwartz's  model  for  biofeedback  (Schwartz,  1979). 


Schwartz's  systems  approach .  Figure  2  shows  the  model  for  biofeedback 
proposed  by  Schwartz  (1979).  Strictly,  this  model  applies  to  biofeedback  of 
peripheral  skeletal  and  autonomic  functions.  However,  it  can  also  be  applied 
to  EEG  biofeedback.  Note  that,  as  in  Mulholland's  model,  an  internal  negative 
feedback  loop  is  hypothesized.  During  biofeedback,  an  external  loop  is  added 
(5),  and  it  is  the  association  between  external  and  internal  feedback  informa¬ 
tion  which  presumably  leads  to  improvement  in  control  of  the  target  behavior. 
However,  it  should  be  noted  that,  since  central  mechanisms  involved  in  internal 
monitoring  and  those  involved  in  processing  the  external  biofeedback  information 
may  not  be  the  same,  such  an  association  may  not  lead  to  enhanced  control  when 
external  feedback  is  disconnected. 

A  new  model  of  biofeedback.  Figure  3  shows  the  model  of  biofeedback  that 
underlies  the  present  research.  It  borrows  from  both  the  Mulholland  and 
Schwartz  models,  but  adds  some  new  concepts  from  Dinnat  (1979)  and  the  authors. 

It  is  proposed  that  a  neuro-generator  (1)  is  responsible  for  the  generation  of 
alpha  activity  (2) .  The  presence  of  alpha  activity  gives  rise  to  observable 
states  (3)  which  can  be  either  provided  externally  (4  -  the  feedback  signal  in 
the  biofeedback  experiment)  or  internally  (5).  Evidence  for  internally  provided 
observable  states  could  come  from  finding  that  subjects  are  able  to  increase 
alpha  aotivity  in  the  absence  of  any  external  feedback.  These  observable  states 
are  monitored  by  a  central  monitoring  process  (6)  which  may  also  be  required  to 
monitor  the  external  environment  (7).  As  a  result  of  monitoring,  the  system 
activates  (8)  alpha  mediators  which  are  connected  to  the  neuro-generator.  The 
activator  is  also  involved  in  providing  resources  for  cognitive  tasks  (9).  Dis¬ 
ruption  of  the  biofeedback  system  can  occur  at  two  levels.  First,  if  environmental 
stimuli  force  the  monitoring  system  to  switch  from  the  observation  of  the  alpha- 


Figure  3-  The  proposed  model  of  biofeedback. 


related  states  then  alpha  production  must  necessarily  suffer.  Second,  if  the 
activating  system  is  involved  in  providing  resources  for  cognitive  tasks,  it 
cannot  at  the  same  time  activate  alpha  mediators.  Under  these  conditions  then, 
the  production  of  alpha  would  suffer.  Furthermore,  the  environmental  stimuli 
may  themselves  produce  a  state  which  is  antithetical  to  alpha  production  (cf. 
Mulholland's  work  on  visual  stimuli  and  the  occipital  EEG). 

Note  that  control  of  alpha  production  occurs  when  the  association  is  made 
between  the  mediators,  the  alpha  state  and  the  external  observable  state  (feedback 

[12] ).  If  the  mediators  are  consciously  generated  (see  below)  and  are  themselves 
observable  internal  states  (5),  then  transfer  of  control  from  the  feedback  to  the 
no-feedback  situition  can  occur  after  internal  and  external  states  are  associated 

(13) -  Controlled  generation  without  feedback  would  then  require  that  an  asso¬ 
ciation  between  alpha  states  and  alpha-related  internal  stimuli  had  been 
developed  (14), 

Mediation 

Sufficient  data  is  now  available  to  indicate  that  at  least  some  individuals 
can  gain  control  of  the  production  of  alpha  activity.  Of  central  concern  here 
is  the  question  of  how  this  control  is  achieved  —  that  is,  what  are  the  mediating 
processes?  In  the  previous  section,  we  noted  that  interference  with  the  pro¬ 
duction  of  alpha  activity  could  be  attributable  to  the  draining  of  resources 
required  for  mediation  by  cognitive  activity.  Thus,  the  extent  of  interference 
would  depend  on  the  degree  to  which  cognitive  and  mediating  activity  draw  on 
the  same  resources. 

The  question  of  mediation  has  been  of  central  concern  to  many  researchers 
in  the  biofeedback  area.  In  general,  two  classes  of  mediation  are  identified 
(Katkin  and  Murray,  1968)  —  somatic  and  cognitive .  Somatic  mediation  involves 
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the  use  of  skeletally  controlled  behavior  to  generate  the  desired  output, 
while  cognitive  mediation  involves  the  use  of  "thoughts." 

The  somatic  meditors  reported  by  subjects  enhancing  alpha  activity  include 
relaxation  (e.g.,  Brown,  1971;  Nowlis  and  Kamiya,  1970)  and  "not-focussing" 

(e.g.,  Nowlis  and  Kamiya,  1970;  Plotkin,  1976).  Relaxation  would  seem  to  be 
especially  susceptible  to  a  variety  of  "noxious"  environmental  events  and 
associated  affective  processes  (such  as  stress  and  anxiety)  which  produce 
increases  in  muscle  tension.  On  the  other  hand,  mediation  involving  "not- 
focussing"  would  be  influenced  by  any  situation  requiring  visual  processing. 

The  cognitive  mediators  include  letting  go,  floating,  awareness  "in  back," 
etc.  (Nowlis  and  Kamiya,  1970).  Although  these  subjective  descriptions  are  vague, 
it  is  not  unreasonable  to  suppose  that  anv  task  or  situation  requiring  cognitive 
activity  would  disrupt  the  cognitive  mediating  process. 

Hypotheses 

*  If  the  revised  model  of  alpha  regulation  was  representative  of  a  true 
volitional  control  system  with  major  cognitive  mediation  components,  then  a 
number  of  hypotheses  could  be  generated  based  upon  this  assumption.  Dinnat 
(1979)  proposed  a  series  of  such  hypotheses,  some  of  which  were  a  direct 
outgrowth  of  his  original  model.  These  were  (1)  that  performance  would  stabilize 
after  sufficient  training  and  that  all  subsequent  performance  would  fall  within 
determinate  bounds;  (2)  in  the  absence  of  feedback  performance  would  eventually 
return  to  baseline  levels;  (3)  performance  with  feedback  would  be  superior  to 
performance  without  feedback;  (4)  if  the  feedback  sensory  mode  inter f erred  with 
alpha  production,  these  effects  would  be  detectable  using  environmental  changes/ 
stimuli  in  that  sensory  mode;  (5)  if  the  feedback  sensory  mode  did  not  interfere 
with  alpha  production,  production  performance  would  be  independent  of  activity 
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in  that  sensory  channel;  (6)  alpha  generation  should  be  subject  to  disruption 
by  enviromental  distractors  in  the  same  way  that  cognitive  tasks  were. 

If  this  last  hypothesis  was  substantiated  one  might  be  able  to  infer  dis¬ 
ruptions  of  cognitive  tasks  by  environmental  stimuli  through  observation  of 
the  disruption  of  alpha  production.  A  series  of  experiments,  each  dependent 
upon  its  immediate  predecessor,  was  conducted  to  evaluate  these  and  other 
hypotheses  concerning  the  production  of  alpha  activity  and  the  extent  to  which 
central  mediation  of  the  response  might  be  involved. 
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PHASE  1 


The  evaluation  of  alpha  generation  as  an  index  of  productivity  potential 
required  development  in  two  major  areas  prior  to  experimentation.  These  two 
aspects  of  system  architecture  were  hardware  devices  and  software  systems . 

Handware  Bsaiasa. 

It  was  desirable  that  the  measurement  system  be  simple,  economical,  and 
transportable  for  possible  future  use  in  the  field.  These  requirements  were 
satisfied  through  the  use  of  a  microprocessor-based  computer  system  with  appro¬ 
priate  input/output  devices  (Figure  4).  This  consisted  of  an  Apple  11+  computer 
with  two  floppy  disk  drives,  12"  c.r.t.  monitor,  and  dot-matrix  impact  printer 
(graphics  capability)  using  a  parallel  interface.  The  system  also  contained 
a  16-channel  analog-to-digital  converter,  RS232  serial  interface,  48K  bytes  of 
memory,  and  a  real-time  clock.  A  smaller  monitor  was  added  to  allow  remote  moni¬ 
toring  of  results  by  experimental  subjects. 

The  microprocessor  system  received  processed  EEG  information  (in  analog 
form)  from  an  Autogen  120  Encephalograph  Analyzer.  It  was  necessary  to  both 
electrically  isolate  this  device  from  the  computer  and  to  scale  down  the  analog 
outputs  of  the  device.  This  was  accomplished  by  use  of  a  simple  and  inexpensive 
optical  interface  constructed  specifically  for  the  project  (see  Appendix  A  for 
schematic).  The  moderate  nonlinearities  (see  response  curves,  Appendix  A) 
inherent  in  such  an  unsophisticated  device  were  easily  compensated  for  in  the 
data  collection  software.  The  two-channel  device  proved  extremely  reliable, 
requiring  only  infrequent  replacement  of  the  batteries  supplying  one  of  the 
two  sources  of  interface  electrical  power. 
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FLOPPY  DISK 
DRIVES 


Figure  M.  Data  collection  system  hardware  configuration. 
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Software  Systems 

All  data  collection  and  manipulation  operations  were  performed  using  the 
microprocessor  system  in  floating-point  BASIC.  Experimental  trials  were  inter¬ 
mediately  stored  in  memory  and  transferred  to  disk  at  the  end  of  each  trial. 

The  joftware  provided  not  only  real-time  graphic  displays  of  performance  but 
also  the  capability  of  graphic  presentation  of  results  at  the  termination  of 
each  trial.  Data  manipulation  software  developed  for  the  project  provided 
hard  copy  of  data  listings,  graphic  results,  and  statistical  analyses.  The 
packages  also  provided  multiple  overlay  capabilities  for  comparative  graphics. 
Software  listings  are  contained  in  Appendix  B. 

Method 

The  first  major  phase  of  experimentation  dealt  with  the  questions  of  per¬ 
formance  stability  as  it  affected  training  requirements  and  the  role  of  feed¬ 
back  in  performance  enhancement.  Ten  subjects,  5  male  and  5  female,  were 
selected  from  among  university  students  responding  to  advertising  for  the  study. 
Although  these  participants  were  guaranteed  minimum  wage  for  the  time  of  their 
service,  most  were  motivated  by  the  potential  of  learning  to  control  their  brain¬ 
wave  activity.  In  this  respect  there  was  some  preselection  of  subjects  inasmuch 
as  they  were  not  conscripted  from  undergraduate  psychology  courses  as  is 
often  the  case.  These  subjects  were  given  a  briefing  on  biofeedback  as  a  general 
concept  and  connected  to  the  measurement  system  using  three  sponge  electrodes. 

The  two  active  sites  measured  were  left  frontal  (above  the  left  eye)  and  mid- 
sagital.  The  reference  electrode  was  placed  behind  the  right  ear  over  the 
mastoid  bone. 

Subjects  were  instructed  to  seat  themselves  comfortably,  relax,  close  their 
eyes,  and  avoid  the  use  of  facial  or  Jaw  muscles  for  a  few  minutes.  After  the 
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subject  appeared  relativ:  .y  confortable  a  two-minute  baseline  trial  was  initiated 
at  a  signal-amplitude  criterion  of  30  microvolts.  At  the  conclusion  of  this 
trial  the  subject  was  instructed  to  open  his/her  eyes  and  look  into  the  center 
of  a  large  piece  of  cardboard  that  was  devoid  of  any  distinguishing  marks.  A 
baseline  was  then  measured  in  this  eyes-open  condition  at  a  criterion  of  30 
microvolts.  This  was  fallowed  by  instruction  on  the  various  strategies  that 
might  be  useful  in  elevating  alpha  production.  Each  subject  then  received  an 
average  of  17  minutes  of  practice  time  with  auditory  feedback  (three  5-minute 
trials  with  interpolated  discussion) .  Feedback  consisted  of  a  tone  varying 
in  frequency  and  amplitude  as  did  the  averaged  EEG  activity  while  within  the 
bounds  of  alpha  (8  to  13  Hz). 

At  the  end  of  this  instructional  period  a  series  of  two-minute  performance 
trials  was  conducted  to  determine  at  what  point  short-term  performance  reached 
an  asymptote.  Criterion  for  reaching  stable  performance  was  set  as  three  con¬ 
secutive  trials  on  which  final  percentage  scores  increased  by  no  more  than  five 
percentage  points.  The  dependent  measure  was  percent  time  in  alpha,  calculated 
as  time  in  alpha  divided  by  total  time  on  task.  This  fraction  was  continuously 
recalculated  throughout  the  trial  at  the  rate  of  six  times  per  second.  The  mean 
number  of  trials  required  to  satisfy  this  criterion  was  3*9  (s.d.s  1.28).  This 
normally  ended  the  first  session  as  overall  time  in  experimentation  was  usually 
one  hour  by  this  time  and  some  signs  of  subject  fatigue  were  often  evident. 

The  next  segment  of  performance  assessment  involved  trials  to  criterion 
with  eyes  open  and  the  effects  of  feedback  withdrawal  on  sustained  performance. 

A  baseline  in  the  eyes-open  condition  was  measured  and  subjects  then  proceeded 
as  in  the  first-session  criterion  trials.  These  measurements  were  followed  by 
two  series  of  four-minute  trials  requiring  alpha  production  without  feedback, 
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both  with  eyes  closed  and  eyes  open.  Subjects  again  used  a  relatively  empty 
visual  field  during  the  eyes-open  trials. 

Rwulta  and.  Discussion 

Figure  5A  shows  seven  superposed  criterion  trials  for  subject  #1.  This  is 
immediately  followed  by  a  point-by-point  plot  of  the  mean  plus  and  minus  one 
standard  deviation  for  these  trials  (Figure  5B).  These  results  are  generally 
representative  of  majority  performance  on  the  task,  with  final  performance 
levels  being  reasonably  close.  There  is,  however,  considerable  variation  at 
the  outset  of  the  trials.  This  is  not  suprising  due  to  the  nature  of  the 
dependent  variable  and  its  statistical  tendency  towards  reduced  variability  as 
time  on  task  increases.  This  artifact  is  illustrated  in  Figure  5C,  showing  a 
plot  of  hypothetical  performance  alternating  from  alpha  to  no-alpha  at  2  Hz 
(fraction  recalculated  at  1  Hz).  The  discontinuity  of  the  function  and  its 
rapid  bounding  are  both  evident,  the  discontinuity  contributing  to  maximized 
variation  in  the  early  seconds  of  the  trial  and  the  limiting  nature  tending 
to  stabilize  the  resultant  in  the  latter  stages.  What  is  interesting  in  the 
performance  data  is  that  although  overall  performance  did  stabilize  with  repeated 
trials,  the  manner  in  which  subjects  achieved  that  final  value  did  not.  The 
final  trial  in  the  series  for  subject  #1  is,  in  fact,  the  top  curve  in  the 
early  seconds  in  Figure  5A.  This  curve  clearly  falls  outside  of  a  standard 
deviation  for  a  good  part  of  the  trail. 

The  fact  that  multiple  paths  exist  leading  to  the  same  overall  performance 
suggests  that  there  is  no  stable  "production  function",  underlining  the  random- 
burst  nature  of  alpha.  It  appears  that  overall  alpha  production  as  measured  over  a 
specified  time  period  is  relatively  constant  under  constant  conditions. 

The  distribution  of  this  activity,  however,  may  vary  across  trials. 


A  second  question  area  involved  the  function  of  feedback  in  facilitating 
performance.  Hypothesis  (2),  the  eventual  decay  of  performance  to  baseline 
levels  after  withdrawal  of  feedback,  was  not  supported.  The  ten  subjects  examined 
were  evenly  divided  between  those  whose  performance  eventually  dropped  below 
baseline  levels  and  those  whose  performance  remained  above  baseline.  There 
had  been  some  concern  that  performance  was  not  returning  to  the  baseline  levels 
because  alpha  mediators  were  being  incorporated  into  each  individual's  normal 
relaxation  behavior.  This  was  examined  for  16  subjects  from  the  Phase  2  section 
of  the  study  by  taking  pre-  and  post-training  baselines  for  each  subject.  A 
one-way  analysis  of  variance  of  percent  alpha  scores  demonstrated  no  reliable 
shift  in  baseline  scores  (F(1,15)  *  3-72,  p  >  .05).  A  similar  analysis  was 
conducted  on  cycle  frequency  into  and  out  of  alpha.  Scores  were  transformed  by 
log(x)  to  better  meet  the  assumptions  of  the  test  (Myers,  1972).  Again,  no  re¬ 
liable  baseline  shift  was  found  (F(1,15)  =  1-37,  p  >  .05).  Thus  training 
did  not  reliably  affect  baseline  unintentional  performance.  (ANOVA  summary 
tables  are  contained  in  Appendix  C.) 

It  had  also  been  hypothesized  that  performance  with  feedback  would  be  superior 
to  performance  without  feedback.  This  question  was  approached  by  comparing  the 
average  scores  from  performance  trials  using  feedback  with  scores  from 
performance  trials  without  feedback.  Only  the  first  two  minutes  of  the  no¬ 
feedback  trials  were  used  so  that  comparisons  would  be  made  over  equal  lengths 
of  trial  time.  (No  reliable  differences  were  found  between  summary  measures  for 
the  first  two  minutes  and  summary  measures  for  the  entire  four  minutes  of  the 
original  trial.) 

The  initial  within-subject  analysis  of  variance  of  percent  alpha  by  both 
feedback  condition  and  eyes  open/closed  indicated  reliable  effects  for  feed- 


back  (F(1,9)  =  6.45,  n.  <  .05),  eyes  open/closed  (F(1,9)  =  21.6,  £  <  .01), 
and  the  interaction  of  the  two  F(1,9)  =  8.33,  Ji  <  .025).  This  was  not 
unexpected  as  it  has  already  been  demonstrated  that  alpha  generation  with 
eyes  open  is  less  than  that  achievable  with  eyes  closed.  The  scores  were 
then  transformed  by  baseline  scaling  to  eliminate  the  effect  of  eyes  open/closed 
and  also  transformed  by  log(x+30)  to  fit  the  linear  model  more  accurately. 

This  analysis  yielded  near-reliable  effects  for  feedback  condition  (F(1,9)  = 
4.73:  F(1,9)  at  p( .05)  =  5.12),  and  the  interaction  of  feedback  condition  with 
eyes  open/closed  (F(1,9)  =  5-03).  The  effect  of  eyes  open/closed,  as  antici¬ 
pated,  was  not  reliable  (F(1,9)  =  1.04). 

The  mean  of  baseline-scaled  scores  in  the  eyes-closed  condition  dropped 
from  6 . 3%  with  feedback  to  2.5>  without  feedback.  In  the  eyes-open  condition 
the  same  respective  means  were  -11. 6^  and  2.6%.  The  means  for  these  condi¬ 
tions  are  depicted  in  Figure  6  (converted  back  to  the  original  units  from  the 
transformed  means) .  This  approach  does  make  interpretation  somewhat  more  diffi¬ 
cult.  More  seriously,  however,  this  result  was  quite  different  from  that 
derived  from  the  unsealed  data  where  reliable  decreases  were  evident  from  eyes- 
closed  to  eyes-open  and  from  feedback  to  no-feedback  conditions.  There  still 
remains  a  reasonable  question  as  to  the  validity  of  the  second-session  base¬ 
lines  taken  prior  to  the  no- feedback  trials. 

Gross  departures  from  normality  and  differing  distributions  across  cells  of 
the  design  made  the  same  analysis  impractical  for  cycle  frequency.  The  mean 
differences,  however,  were  only  .17  Hz  at  the  greatest  for  the  raw  data  and 
.01  Hz  for  the  baseline-scaled  data.  There  did  not,  as  such,  appear  to  be  any 
consistent  or  meaningful  differences  in  cycle  frequency  as  a  result  of  feedback 
or  eye-condition  manipulations.  It  should  be  noted  that  no  differences  attri- 
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butable  to  sex  were  detected  for  the  10  subjects  in  Phase  1. 

One  possible  explanation  for  the  eyes-open  increase  from  feedback  to  no¬ 
feedback  conditions  relies  on  feedback  mode  as  a  basis.  Given  the  proposed 
model,  positive  feedback  vould  require  the  cognitive  monitor  to  switch  from 
monitoring  internal  medi-^ors  to  processing  external  feedback.  The  resultant 
lapse  in  mediator  supervision  could  then  result  in  decrements  in  alpha  perfor¬ 
mance.  It  is  not  clear  how  this  process  might  have  have  an  inverse  effect  when 
the  eyes  are  closed.  Phase  2  experimentation  was  directed  towards  further 
examining  this  problem. 


>p  -  5 
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Figure  6.  Retransformed  condition  means:  Baseline-scaled  percent  alpha  by  feed¬ 
back  and  eye  conditions. 
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PHASE  2 


The  second  phase  of  experimentation  was  designed  to  investigate  the  effects 
of  feedback  mode  (auditory/visual)  and  feedback  polarity  (positive/negative) 
upon  alpha  production.  It  has  already  been  mentioned  that  feedback  provided 
through  a  non- inter f erring  sensory  channel  would  more  readily  enhance  performance. 
Visual  activity  is  known  to  inhibit  alpha  activity,  suggesting  that  auditory 
or  tactile  channels  are  better  for  the  transmission  of  feedback  information. 

Most  tasks  of  interest  where  this  metric  would  be  applied,  however,  have  some 
visual  components  and,  as  demonstrated  earlier,  there  may  exist  some  interaction 
between  visual/auditory  activity  and  performance. 

Two  scenarios  can  be  developed  using  the  proposed  model  operating  under 
either  positive  or  negative  feedback  respectively.  Under  positive  feedback 
(stimulus  on  when  alpha  produced)  the  cognitive  monitor  would  be  called  upon 
to  monitor  the  feedback  stimulus  when  it  appeared,  drawing  the  monitor  away 
from  the  mediator  information.  In  the  absence  of  this  input,  the  generating 
processes  could  wander,  leaving  the  criterion  region  of  performance.  The 
external  feedback  would  then  cease,  allowing  the  monitor  to  again  receive 
mediator  information,  completing  the  cycle.  This  behavior  is  in  accord  with 
what  Mulholland  (1977)  had  shown  in  his  report  of  reduced  cycling  variation. 

Negative  feedback  (stimulus  on  when  no  alpha  produced)  may,  however,  have 
a  different  function.  If  we  begin  in  a  no-alpha  state,  the  stimulus  is  on, 
again  drawing  the  monitor  to  external  events  and  away  from  mediator  Information. 
The  only  way  to  enter  criterion  alpha  production  would  be  through  the  normal 
baseline  cycling  from  no-alpha  to  alpha,  turning  off  the  external  stimulus. 


Once  this  was  accomplished  the  monitor  could  then  return  to  receiving  mediator 
information  and  could  remain  in  this  state  for  as  long  as  it  was  not  called  upon 
to  receive  other  information.  This  could  result  in  higher  variability  of  cycling. 
Its  projected  effects  upon  actual  time  in  alpha  and  cycle  frequency  are  unclear 
and  were  an  object  of  investigation. 

An  additional  19  subjects  who  responded  to  the  solicitation  for  participants 
were  screened  for  intentional  alpha  production.  The  procedures  were  the  same  as 
those  used  in  the  initial  stages  of  Phase  1;  initial  baselines,  training, 
secondary  baselines,  and  eyes-closed  intentional  production  with  feedback. 

Six  of  these  individuals  demonstrated  some  ability  to  enhance  their  alpha 
production.  These  individuals  returned  for  a  second  session  during  which 
additional  baselines  were  measured  first.  Two-minute  performance  trials  were 
then  conducted,  in  randomized  orders,  in  each  of  the  conditions  delineated  by 
a  two-by-two  factorial  within-subject  design.  Each  initial  change  in  feedback 
modality  was  accompanied  by  2  minutes  of  practice  in  that  mode  prior  to  perfor¬ 
mance  assessment.  The  independent  variables  were  feedback  modality  (auditory/ 
visual)  and  feedback  polarity  (positive/negative).  The  dependent  variables  were 
percent  time  alpha  produced  and  cycle  frequency. 

A  modification  of  the  original  analog-to-analog  interface  provided  for  either 
visual  (L.E.D.’s)  or  auditory  (digitally  generated  tone)  feedback  (positive  or 
negative)  (schematics  in  Appendix  A).  Auditory  feedback  was  initially  presented, 
in  this  phase,  through  light  earphones  that  fitted  inside  the  outer  ear  (in  Phase 
1  this  was  accomplished  using  a  speaker  near  the  subject) .  These  were  later 
replaced  by  acoustic-tube  earphones  so  that  the  electromagnetic  driving  element 
would  not  be  in  proximity  to  any  of  the  electrodes.  Visual  feedback  was  presented 
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by  use  of  two  red  L.E.D.’s.  These  were  mounted  on  the  inner  surface  of  a  trans¬ 
lucent  industrial  face  shield  and  adj tsted  to  eye  level.  This  provided  a  homo¬ 
geneously  illuminated  background  upon  which  the  stimulus  could  be  detected  and 
minimized  opportunities  for  visual  search  activity. 

Results  and  Discussion 

Two-way  analyses  of  variance  were  conducted  on  raw  and  baseline-scaled 
percent  time  in  alpha  (J  alpha)  and  cycle  frequency  by  feedback  mode  and  feedback 
polarity.  All  untransformed-data  analyses  exhibited  F-ratios  considerably  less 
than  1.  Analyses  performed  on  transformed  scores  (Appendix  C)  still,  at  the  best, 
exhibited  similar  results.  Closer  examination  of  the  within-cell  distributions 
indicated  that  normality  assumptions  had  been  violated.  Transformations  did  not 
help  because  there  was  serious  heterogeneity  among  the  distributions.  Had  the 
distributions  been  nonnormal  but  similar,  distortions  of  the  test  would  have 
been  minimized. 

In  view  of  this  problem  the  Friedman  two-way  analysis  of  variance  by  ranks 
(Horowitz,  1974)  was  used  to  examine  the  data.  The  sample  was  large  enough  that 
the  sampling  distribution  of  the  Friedman  test  followed  the  chi-square  distri¬ 
bution  with  3  degrees  of  freedom.  The  critical  value  for  p(.05)  was  7.81. 

Obtained  values  were  2.6  for  %  alpha  and  2.4  for  cycle  frequency.  The  baseline- 
scaled  computations  required  special  chi-square  (r)  tables  for  a  smaller  n 
(Siegel,  1956),  giving  a  critical  value  for  p'-05)  of  7.5.  Obtained  values 
were  0.9  for  baseline-scaled  %  alpha  and  3-9  for  baseline-scaled  cycle  frequency. 
Thus,  no  reliable  effects  due  to  treatments  were  found.  This  is  undoubtedly  due 
to  the  restricted  sample  size  and  high  individual  variability. 

Some  other  apparently  consistent  relationships  had  been  observed  in  the 
summary  data  and  a  decision  was  made  to  investigate  the  correlation  between  the 
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dependent  variables.  The  overall  correlation  (Appendix  D)  across  all  subjects 
and  all  conditions  was  R  =  .71  [t(406)  s  20.35;  p  <  .001].  This  plot  exhibited 
a  marked  hook,  however,  the  curve  doubling  back  to  high  %  alpha,  low  frequency. 
Subsequent  separation  of  the  data  by  eyes-closed  (1,3, 5, 7, 8)  and  eyes-open 
(2,4,6,9110)  conditions  showed  that  these  nonlinear  effects  were  concentrated 
in  the  eyes-closed  trials  [R  s  .596;  t(245)  =  11.63*  P  <  .001].  The  eyes- 
open  trials  exhibited  a  very  strong  linear  relationship  between  the  variables 
[R  s  .88;  t ( 159)  =  23*09,  p  <  .001]  (  see  Appendix  D) .  Although  these 
correlations  were  also  examined  by  each  of  ten  conditions,  microvolt  levels, 
and  subjects,  eyes-open  versus  eyes-closed  appeared  to  be  the  greatest  defining 
contrast.  This  was  undeniably  due  to  the  alpha  blocking  found  with  visual 
activity,  reducing  variability  in  that  condition. 

This  correlation  between  the  dependent  variables  has  implications  for 
multivariate  metrics  that  were  earlier  thought  possible.  Preliminary  data 
collected  on  and  analyzed  from  strip-chart  recordings  suggested  that  these 
two  variables  had  a  certain  measure  of  independence.  Examination  of 
activity  during  various  tasks  (mental  arithmetic,  reading,  visual  tracing, 
monitoring  spoken  text  for  target  words)  suggested  that  some  tasks  that  could 
not  be  differentiated  by  %  alpha  alone  could  be  differentiated  by  concurrent 
examination  of  cycle  frequency.  If  the  high  positive  correlation  of  these 
two  variables  extends  beyond  the  confines  of  this  particular  experimental 
paradigm,  then  the  use  of  these  variables  in  a  multivariate  discrimination  scheme 
may  not  be  realized. 
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CONCLUSIONS 

Although  the  model  of  biofeedback  and  electrocortical  activity  control 
proposed  herein  was  consistent  with  previous  research  findings,  additional  con¬ 
firmation  of  its  validity  was  not  obtained  from  the  present  study.  It  is  diffi¬ 
cult  to  determine  which  factors  were  most  influential  in  producing  this  outcome. 
Certainly,  the  presence  of  tremendous  variability  between  and  within  subjects 
coupled  with  the  low  incidence  of  alpha  production  enhancement  makes  this  par¬ 
ticular  metric  an  unattractive  one.  While  it  appears  to  have  a  high  degree  of 
sensitivity  to  enviromental  changes,  it  is,  perhaps,  overly  sensitive  and  thus 
unstable.  More  productive  metrics  for  the  evaluation  of  environmental  impacts 
on  behavior  might  be  developed  from  evoked  cortical  potentials,  physiological 
measures,  or  some  combination  thereof. 
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APPENDIX  A: 
SCHEMATIC  DIAGRAMS 
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Figure  7.  Schematic  diagram:  Analog-to-analog  interface,  channel  1. 
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ANALOG-TO-ANALOG  OPTICALLY- ISOLATED  INTER¬ 
FACE:  AUTOGEN  I20A  ELECTROENCEPHALOGRAPH 
OUTPUT  (GROUND-REFERENCED  CHANNEL) TO 
APPLE n  COMPUTER  A“TO-D  BOARD. 
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Figure  8.  Schematic  diagram:  Analog-to-analog  interface,  channel  2. 
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Figure  10.  Schematic  diagram:  Digital  stimulus  controller. 
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Figure  11.  Schematic  diagram:  Digital  tone  generator. 
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APPENDIX  B: 
SOFTWARE  LISTINGS 


3PRINT  DATA  COLLECTION  AND  MANIPULATION  PROGRAM 

7SYNTAX  ERROR 
3LIST 

I  DIM  PC2(  300  ),SUZ(  300  )»FCOUNTZ(  300  60  ) 

3  DSTART  =  24586 5USTART  =  25500: A  =  -  158721  REM  A  VALUE  FOR  A/D  IN  SLOT 

n 

4  PRINT  :  PRINT  "LOADING  WRITE  FLAGS" 

5  HIMEM:  16383 ID*  =  CHR$  <4)5SSTART  =  25572 

6  PRINT  D$»"BLCAD  UFLAG, A25500" :  REM  SOFTWARE  PROTECTION  TO  PREVENT  OVER 

WRITES 

7  SBNUM  =  0 

8  CALL  62450:  TEXT  5  HOME  IX  =  FRE  (0) 

9  PRINT  "1  FC:  DATA  COLLECTION." 

10  INPUT  “2  FOR  DISK  RETRIEVE.  "?G 

II  IF  G  =  2  THEN  GOTO  3000 

12  GCSUB  9000 

13  HGR2 

14  SAMPS  =  CJDEXER  =  OISCREENS  =  OISY  =  OJIA  =  OJSl  =  OJCOUNT  =  OIOFLAG  = 

15  CALL  62450 

16  XM  =  I!  GCSUF  7000 

18  FOR  I  =  10  TO  270 

19  SFLAG  =  0: SAMPS  =  SAMPS  +  1 

20  POKE  A  +  1,6 

21  IN  =  PEEK  (A)I  IF  IN  <  85  AND  IA  =  0  THEN  GOTO  20 
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22  POKE  A  +  1.7:  IF  IA  =  0  THEN  IA  =  1 


23  PER  =  PE 

24  GOSUB  10 
*»«?  pv  =  no 


PEEK  (A) 


25  PY  =  128  -  (PER  /  2) 

26  OY  =  128  +  (( IN  -  75)  /  50) 

27  IF  IN  >  85  THEN  SI  =  SI  +  1ISFLAG  =  1 

28  SY  =  128  -  ( 128  *  (SI  /  (( I  -  9)  +  (SCREENS  *  261)))) 

29  IF  SY  >  191  OR  SY  <  0  THEN  GOTO  300 

30  HCCLOR=  7*.  HPLOT  I,PY 

31  HPLOT  1,123  TO  ItOY 

32  HPLOT  I,SY 

33  IF  SFLAG  =  OFLAG  THEN  GOTO  40 

34  COUNT  =  COUNT  t  II OFLAG  =  SFLAG 
40  IF  SAMPS  =  6  THEN  GOSUB  190 

50  NEXT  1 1 SCREENS  =  SCREENS  +  1 

51  IF  SCREENS  <  3  *  MINS  THEN  GOTO  15 

52  GOTO  1000 

100  IF  PER  <  28  THEN  PER  =  PER  -  185  GOTO  110 

101  IF  PER  <  193  THEN  PER  =  PER  -  131  GOTO  110 

102  IF  PER  <  213  THEN  PER  =  PER  -  10J  GOTO  110 

103  IF  PER  <  221  THEN  PER  =  PER  -  5:  GOTO  110 

104  IF  PER  <  230  THEN  GOTO  110 

105  IF  PER  <  237  THEN  PER  =  PER  +  5:  GOTO  110 

106  PER  =  PER  +  10 

107  REM  SOFTWARE  COMPENSATION  FOR  NONLINEAR  ANALOG  INPUT 
110  RETURN 

190  D1  =  INT  (PER  /  2.55 ):D2  =  (PER  /  2.55) 

195  IF  (02  -  01)  >  .49  THEN  D1  =  D1  +  1 

200  PCX(DEXER)  =  D1 

203  D 1  =  INT  ((128  -  SY)  /  1.28)102  -  (( 128  -  SY )  /  1.28) 

206  IF  (D2  -  01)  >  .49  THEN  D1  =  D1  +  1 

210  SUZ(DEXER)  =  B1 

220  IF  DEXER  =  0  THEN  FCOUNTZ(  DEXER )  =  COUNT.*  GOTO  230 
225  FCOUHTft DEXER )  =  COUNT  -  OCOUNT 

230  DEXER  =  DEXER  +  1 

231  OCOUNT  =  COUNT 
235  SAMPLE  =  0 

240  IF  DEXER  >  (60  *  MINS)  -  1  THEN  GOTO  1000 
250  RETURN 
300  TEXT 


PRINT  "SY  =  "iSY 
PRINT  "SCREENS  =  "JSCREENS 
INPUT  "2  TO  CONTINUE,  1  TO  QUIT"iG 
IF  G  =  2  THEN  GOTO  13 
CALL  62450*.  TEXT 
PRINT  CHR*  (7):  PRINT  CHR*  (7) 

PRINT  .*  PRINT  "TRIAL  CONCLUDED":  PRINT 

PRINT  "(1)  LIST  DATA  (TABLE/GRAPH)":  PRINT 

PRINT  "(2)  RECOLLECT  DATA  (NO  STORAGE  )“J  PRINT 

INPUT  "(3)  SEND  DATA  TO  DISK  “JG 

IF  G  =  2  THEN  CALL  62450:RC  =  1J  GOTO  12 

IF  G  *  3  THEN  GOTO  2000 

GOTO  3092 
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m 
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2052 
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2062 

r>r>Ai 


CALL  624501  TEXT 
GOSUB  9800 

INPUT  "LIST(L)  OR  NEXT  IB  <N)?”.*A* 

IF  AS  =  “N“  THEN  Z  =  Z  +  11  GOTO  3013 
BZ  =  0 

CZ  =  FZ  4  19 
FOR  I'EXER  =  BZ  TO  CZ 
PRINT  PCX(  DEXER  ) * “  "»SUX< BEXER).” 

NEXT  DEXER J  PRINT 

INF'u'  “RETURN  TO  CONTINUE” )A*:  PRINT 
+  20 

60  *  MINS  THEN  GOTO  1059 
PRINT-LISTING  ENDED”:  PRINT 
1  TO  LIST  AGAIN” 


.FCOUNTZ( BEXER) 


BZ  = 

IF  i. 

prin; 

PRINT 


PRINT  “2  TO  RUN  PROGRAM  AGAIN” 


PRINT 
INPUT 
IF  G  = 

IF  G  = 

IF  G  = 

GOTO  8 
CALL  624 
GOSUB  9800 
PRINT 
HOME 


TO  SEND  TO  DISK" 

TO  CLEAR  RECORD  PROTECTION” iG 
THEN  GOTO  3092 
THEN  GOTO  2000 
THEN  GOSUB  5000*.  GOTO  1150 


INPUT  “RETURN  TO  SEND  TO  DISK'.AS 

_  PRINT  “CHECKING  FILES...” 

FOR  I  =  0  TO  5? 

UF%< I )  =  PEEK  ( USTART  +  I ) 

NEXT  I 
I  =  0 

IF  UF7.(  I )  =  0  THEN  GOTO  2015' 

I  =  I  +  lJ  IF  I  >  59  THEN  GOTO  2011 
GOTO  2003 

PRINT  “ALL  RECORDS  FILLED” 

INPUT  A**.  GOTO  1000 

wfzci)  =  i:z  =  i:  poke  (wstart  +  n»i 

CALL  62450:  HOME  :  PRINT  "SENDING  TO  FILE  B 
FOR  I  =  0  TO  (60  *  MINS)  -  1 
+  (3  *  I))»PCZ(I) 


I))»SUZ(  I) 

m.FcouNTzm 


POKE  ( DSTART 
POKE  (DSTART  t  1  M3  * 

POKE  (DSTART  +  2  +  <3  * 

NEXT  I 
POKE  24576 .SBNUM 
POKE  24577 ,CDN 
POKE  24573. TL 
POKE  24579. DAY 
POKE  24580. MTH 
POKE  24531. YEAR 
POKE  24532. HOUR 
POKE  24583. MNTS 
POKE  24534 .(MINS  *  10) 

POKE  24535. M VOLTS 

PRINT  D$;"BSAVE  B" ;Z» ” .A” .DSTART  -  10.” »L” {< 180  *  MINS)  + 
PRINT  B$» ’BSAVE  WFLAG.A25500.L60* 

PRINT  :  PRINT  "TRANSMISSION  COMPLETE":  PRINT 
PRINT  "1  TO  RETRIEVE  AND  LIST." 

INPUT  "2  TO  COLLECT  DATA  AGAIN. "5G 
IF  G  =  2  THEN  GOTO  8 
CALL  62450 
PRINT 

PRINT  *.  PRINT  ”1  FOR  SINGLE  RECORD” 

PRINT  ”2  FOR  OVERLAYS  " 

INPUT  ”3  FOR  STAT  OPTIONS  ”iG 
IF  G  =  2  THEN  GOTO  6000 

IF  G  =  3  THEN  GOTO  3500 

INPUT  “FILE  NUMBER  <0-5?){  ”.Z 
PRINT  :  PRINT  "ACCESSING  FILE  B'.Z 
GOSUB  9500 

IF,ER  *  1  THEN  PRINT  "FILE  EMPTY" JER  =  O’.  GOTO  3012 

tel  mNT 

INPUT  ”2  FOR  GRAPH." 5G 
IF  G  *  1  THEN  GOTO  1050 

HGR2  :  HCOLORa  7 

XM  =  0:  GOSUB  7000 
G  =  o:  GOSUB  8000 
INPUT  as:  CALL  62450J  HOME  *.  TEXT 
INPUT  "PRESS  RETURN  FOR  CYCLE  FREQUENCY  -fAS 

X«Vo:  GOSUB  7000 
FCR  I  *  0  TO  (60  *  MINS)  -  1 


10 
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3560 

3565 

3570 

3575 

3580 

m 

3595 

3600 

3605 

3610 

3615 

3620 

3621 

3622 
3625 
3630 
3635 
3640 
3645 
4000 
4002 
4005 
4010 
4015 
4020 
4025 
4030 
4035 
4040 
4045 
4047 
4050 
4055 
4060 

4065 

4070 

4075 

4080 

4085 

4090 

4095 

4100 

4105 

4110 

4115 

4120 

4125 

4130 

4135 

5000 


51) 

52) 


FR  =  FCCUNT%<  I ) 

HC01QR=  21  HP10T  I  4  10,128  TO  I  4  10,128  -  (6,4  *  FR) 
NEXT  I 

INPUT  Ail  HOME  :  TEXT 
CALL  62450:  GOTO  1150 
REM  *STAT  ROUTINE* 

HOME  J  PRINT  “PROCESS  FILES  OPTION:' 

PRINT  :  PRINT  “(1)  SHOW  TRIAL  AND  MEAN,  SD  * 

PRINT  :  INPUT  “(2)  ADD  FILES  TO  STAT  BASE’ )G 
IF  G  =  2  THEN  GOSUB  4000:  GOTO  3004 
HGR2  :  HCOLOR=  7 IX M  =  0:  GOSUB  7000 
HCOLOR=  i:  GOSUB  8000:  HCOLQR=  2 
PRINT  DiJ’BLOAD  BSTAT,A25361" 

FOR  I  =  0  TO  (60  *  MINS)  -  1 

N  =  PEEK  ( SSTART  4(3*  I)>:  IF  N  =  0  THEN  GOTO  3610 
XSUM  =  PEEK  (SSTART  4  1  4  (3  *  I)) 

X2SUII  =  PEEK  ( SSTART  4  2  4(3*1)) 

XBAR  =  (XSUM  /  N>  *  10 

SDEV  =  10  *  <  SQR  ((N  *  X2SUM )  -  (XSUM  t  2)))  /  N 

51  =  XBAR  4  SDEV.  IF  SI  >  100  THEN  SI  =  100 

52  =  XBAR  -  SDEU:  IF  S2  <  0  THEN  S2  =  0 

HPLQT  I  4  10,128  -  <1.28  * 

HPLOT  I  4  10,128  -  (1.28  * 

NEXT  I 

INPUT  A*JX  =  FRE  (0) 

CALL  624505  HOME  :  TEXT 

PRINT  M  >  TO  ADD  THIS  TRIAL  TO  SUMMARY  STATS" 

’  (LAST  FILE  ADDED  =  ")  PEEK  (25561))“  )“ 

’  (THIS  FILE  IS  riZ,"  )" 
t  PRINT  “(2)  TO  RETURN  TO  LISTING  OPTIONS" 

INPUT  ’(3)  TO  RETURN  TO  DATA  COLLECTION  “,*G 
6  THEN  GOTO  8 

IF  G  =  2  THEN  GOTO  3004 

RC  =  i:  GOSUB  4000:  GOTO  8 
REM  *ADD  FILES  TO  STAT  DATABASE* 

IF  RC  =  1  THEN  RC  =  OJTR  =  Z:UR  =  Z:  GOTO  4035 
PRINT  :  PRINT  'ADD  TO  STAT  DATABASE" 

PRINT  :  INPUT  "1ST  FILE:  “ )TR 
PRINT  J  INPUT  “LAST  FILE  1  ‘JUR 
IF  UR  <  TR  THEN  UR  =  TR 
PRINT  “ACCESSING  DATA  BASE’ 

PRINT  I<$>  "BLOAD  BSTAT»A25561" 

FOR  K  =  TR  TO  UR 

PRINT  "ACCESSING  FILE  "iKU  =  KJ  GOSUB  9500 
IF  ER  -  1  THEN  PRINT  "EMPTY  FILE  ENCOUNTERED" 

PRINT  “FILE  " >K» "  LOADED" 

FOR  J  =  0  TO  (60  *  MINS)  -  1 


PR  INI 
PRINT  " 
PRINT  : 
PRINT  : 
IF  G  = 


PRINT  :  GOTO 


<3  *  -J>),<  PEEK  (SSTART  4  (3  *  J>>)  +  l 
1  4  (3  *  J)),(  PEEK  (SSTART  4  1  4  (3  t  J)))  4 


6020 

6025 

6040 


POKE  (SSTART  4 
POKE  (SSTART  4 
/  10) 

POKE  (SSTART  4  2  4  (3  *  J)M  PEEK  (SSTART  4  2  4  (3  *  J)))  4 
)  /  10)  t  2) 

NEXT  Jl  PRINT  ’FILE  ',*K»'  ADDED’ J  PRINT 
NEXT  K 

PRINT  "SAVING  LAST  I.D." 

POKE  2556 1,K 

FOR  I  =  24576  TO  24585 

POKE  (I  4  986),  PEEK  (I) 

NEXT  i:  PRINT 

PRINT  "RETURNING  TO  STAT  FILE’ 

PRINT  D*)"BSAVE  BSTAT,A25561 ,L911’ 

PRINT  .*  PRINT  ’COMPLETED:’ 

PRINT  ’(1)  TO  ADD  MORE’ 

INPUT  "(2)  TO  EXIT  ",*G 
IF  G  =  1  THEN  GOTO  4005 
RETURN 

REM  CLEAR  FILE  PROTECTION 

PRINT  .*  PRINT  ‘CLEAR  FILE  PROTECTION:* 

PRINT  :  INPUT  “FROM  (1ST  FILE):  “)G1 
PRINT  t  INPUT  'TO  (2ND  FILE):  "JG2 
IF  G2  <  Gl  THEN  I  =  Gil Gl  =  G2JG2  =  I 
FOR  I  =  Gl  TO  G2 
POKE  (USTART  4  I),0 

NEXT  i:  PRINT  DiJ’BSAVE  UFLAG»A25500,L60" 

CALL  62450:  HOME 
RETURN 

■  '  -  variables: 


PRINT 

PRINT 

PRINT 

PRINT 


PRINT 

PRINT 

INPUT 

INPUT 


1  =  OVERALL  X" 

2  =  INSTANTANEOUS  X" 
’  3  =  CYCLE  FREQUENCY 
’1ST  FILE.'  NT R 


1  »G 


4080 

<SUZ<J) 
( ( SUX( J 


PRINT  :  INPUT  -2ND  FILE:  "*UR 
PRINT  :  PRINT  'ACCESSING  VARIABLE  "*G 
PRINT  “  FILES  “*TR*“  4  "*UR 
Z  =  TR 
GOSUB  9500 

IF  ER  =  1  THEN  PRINT  "FILE  EMPTY" JER  =  0!  GOTO  6040 
HGR2  :  CALL  62450 
HC0LOR=  7 5XM  =  0:  GOSUB  7000 
HCOLQR=  It  GOSUB  8000 
Z  =  UR:  GOSUB  9500 

IF  ER  =  1  THEN  TEXT  5  PRINT  "FILE  EMPTY" JER  =  05  GOTO  6040 
HCOLQR=  2!  GOSUB  8000 
INPUT  A$:  CALL  62450:  HOME  :  TEXT 
PRINT  "1  FOR  MORE  OVERLAYS," 

INPUT  "2  TO  EXIT.  '5G 
IF  G  =  1  THEN  GOTO  6000 
GOTO  8 
HCOLOR=  7 

HPLOT  5*0  TO  5f 131  TO  270*131 
FOR  I  =  0  TO  10 

J  =  35  IF  I  =  0  OR  I  =  5  OR  I  =  10  THEN  J  =  0 
HPLOT  J*12.8  *  I  TO  5*12.8  *  I 
NEXT  I 

IF  XM  =  1  THEN  GOTO  7070 
FOR  I  =  0  TO  3 
FOR  J  =  0  TO  5JK  =  2 
IF  J  =  0  THEN  K  =  7 
IF  J  =  3  THEN  K  =  4 
QX  =  (I  *  60)  +  (J  *  10)  +  10 
HPLOT  QX *131  TO  QX*131  +  K 
NEXT  J 
NEXT  I 
RETURN 

FOR  I  =  0  TO  (60  *  MINS)  -  1 
IF  G  =  1  THEN  GOTO  8010 

IF  G  =  2  THEN  GOTO  8020 

IF  G  =  3  THEN  GOTO  8030 

IF  G  =  0  THEN  HCCLOR=  1 

SY  =  128  -  (  1.28  *  SUX(I)) 

IF  SY  <  0  THEN  SY  =  0 
HPLOT  I  +  10*SY 
IF  G  >  0  THEN  GOTO  8050 

IF  G  =  0  THEN  HCOLQR=  2 

PY  =  128  -  (1.28  *  PCX(I)) 

IF  PY  <  0  THEN  PY  =  0 
HPLOT  I  +  10*PY 
IF  G  >  0  THEN  GOTO  8050 

IF  G  =  0  THEN  HCOLOR=  7 

FR  =  FCOUNT%(  I ) 

IF  G  =  0  THEN  HPLOT  I  +  10,128  TO  I  +  10,128  -  (FR){  GOTO  8050 
HPLOT  I  +  10,128  -  (6.4  *  FR) 

G  =  G 
NEXT  I 
RETURN 

REM  SUBROUTINE  FOR  ID  &  TIME 
IF  RC  =  1  THEN  RC  =  0:  GOTO  9050 
PRINT  5  INPUT  "SUBJECT  NUMBER*.  "*SBNUM 
PRINT  :  INPUT  “CONDITION:  “ *CBN 
PRINT  5  INPUT  'TRIAL:  "*TL 
PRINT  :  INPUT  'TRIAL  LENGTH  (MINUTES):  "*MINS 
IF  MINS  >  4  THEN  MINS  4 
IF  MINS  <  1  THEN  MINS  =  1 
PRINT  :  INPUT  "AMPLITUDE  IN  MVOLTSJ  'JMVOLTS 
PRINT  :  INPUT  "HIT  RETURN  TO  START" *A$ 

REM  GET  THE  TIME  &  BATE 
PRINT  D$,"IN*4"J  REM  INPUT  TO  CLOCK 
INPUT  '  "JT$:  REM  GET  TIME 
PRINT  D<5nIN*0":  REM  IN  TO  KEYBOARD 
MTH$  =  LEFT!  (T$*2) 

BAYS  =  MID*  ( Tt*4*2  ) 

HOURS  =  MIDS  (  T5*7 *2  ) 

MNTSS  =  MIPS  (  Tt* 10*2  ) 

MTH  =  VAL-(MTHS):  REM  DEC  FROM  STRING 
DAY  =  VAL  (-DAYS) 

HOUR  =  VAL  (HOURS) 

MNTS  =  VAL  ( MNTSS  ): YEAR  =  80 
RETURN 

REM  SUBRT  TO  RETRIEVE  A  FILE 
REM  NEEDS  “Z"  FROM  MAIN  PROGRAM 
IF  PEEK  (USTART  +  Z)  =  0  THEN  Eft  *  1 S  GOTO  9580 


DEC  FROM  STRING 
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« 

d 

I 

1 

5 

1 


m 

9512 

9513 

9514 

9515 

9516 

9517 

9518 

9519 

9520 
9530 
9540 
9550 
9560 
9570 
958 C 
9800 
9805 
9810 

m 

9833 

9836 

9840 

9850 

9360 

9866 

9868 

9870 

9880 


PRINT  IWBLDAD  B'SZJ* 

SBNUM  =  PEEK  (24576) 

CBN  s  PEEK  (24577) 

TL  =  PEEK  (24578) 

DAY  =  PEEK  (24579) 

NTH  =  PEEK  (24530) 

YEAR  =  PEEK  (24591) 

HOUR  =  PEEK  (24582) 

MNTS  =  PEEK  (24583) 

KINS  =  (  PEEK  ( 24584 ) )  /  10 
KVOLTS  =  PEEK  (24585) 

FOR  I  =  0  TO  (60  *  MINS)  -  1 
PC%(  I )  =  PEEK  ( BSTART  +  (3  * 


> A* i BSTART  -  10 


I)) 


(3  t  I>> 

’  +  (3*1)) 


SU%(I)=  PEEK  (BSTART  +  1  4 
FCOUNTX(I)  =  PEEK  ( DSTART  4 
NEXT  I 
RETURN 

REM  IB  DISPLAY  SUBROUTINE 
PRINT  "RECORD*  "?Z 
PRINT  "SUBJECT  *  “ 5 SBNUM 
PRINT  “DATE!  "iMTH»M/*»BAY?V"tYEAR 


PRINT 
PRINT  .* 
PRINT  : 
PRINT 
IF  MNTS 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT  AS 
RETURN 


>  9  THEN  PRINT  "TIME:  ■ 
TIKEt  "♦ HOUR? "*0“i MNTS 
PRINT  "CONDITION:  “JCBN 
TRIAL  NUMBER.’  “JTL 


IHOURiTJMNTS!  GOTO  9840 


LENGTH:  " JMINSi "  MINUTES* 

LEVEL  =  " 5MVOLTS5"  MICROVOLTS* 


I 

I 

i 


2PRQGRAM  TO  GENERATE  HARD  COPY  OF  DATA  GRAPHICS 


’SYNTAX  ERROR 
jLIST 

1  Hi  =  CHRi  (4) 

2  CALL  62450:  PRINT  "LOADING  FILE  INFORMATION" 

3  DSTART  =  245865USTART  =  25500ISSTART  =  25573 

5  HIwEik*#  16384 

6  PRINT*  *Dft»"BLOAD  UFLAGf  A25500»D2"  5SBNUM  =  0 

8  CALL  62450:  TEXT  HOME  JX  =  FRE  (0)5  GOTO  3000 

1050  CALL  62450:  TEXT  !  GOSUB  9800 

1051  PRINT  5  PRINT  "SAMPLE  /  ICST%  /  CUMZ  /  FREQ  PRINT 
1055  BZ  =  0 

1059  CZ  =  BZ  +  19 

1060  FOR  BEXER  =  BZ  TO  CZ 

1061  A1  =  PEEK  (DSTART  +  (3  ft  DEXER >)5A2  =  PEEK  ([START  +  60  +  ( 3  ft  DEXE 
R )  > 

1062  Bl  =  PEEK  (DSTART  +  1  +  (3  ft  DEXER)):B2  =  PEEK  (DSTART  r  61  +  <3  ft 
DEXER ) ) 

1063  Cl  =  PEEK  (DSTART  +  2  +  <3  ft  DEXER ))5C2  =  PEEK  (DSTART  +  62  +  < 3  ft 
DEXER ) ) 

1065  IF  CP  =  0  THEN  GOTO  1072 

1070  PRINT  SPC(  2);BEXER  +  l?")“f  SPC(  5)fAli  SPC(  6)»B1?  SPC(  6)fCli  SPC( 
12  )* DEXER  +  615"  )"  5  SPC(  5)5A25  SPC(  6)?B25  SPC(  6)5C2 

1071  GOTO  1080 

1072  PRINT  3PC(  2 ) 5 DEXER  +  1>")U*  SPC(  5)5A15  SPC(  6)5Bli  SPC(  6  )»C1 

1080  NEXT  DEXER 

1081  IF  CP  =  1  THEN  GOTO  1091 

1082  PRINT 

109.0  INPUT  “RETURN  TO  CONTINUE"  5Ai 
1091  BZ  =  BZ  +  20:  IF  CP  =  0  THEN  PRINT 

1095  IF  CP  =  1  THEN  GOTO  1120 

1096  IF  BZ  <  60  *  MINS  THEN  GOTO  1059 

1097  GOTO  1140 

1120  IF  BZ  *  (60  *  (MINS  -  1))  THEN  GOTO  1140 
1122  IF  BZ  =  60  THEN  BZ  =  120 

1124  IF  BZ  =  120  THEN  PRINT  CHRi  (12):  PRINT  :  PRINT  l  PRINT  5  PRINT  SPC< 
70 )» "PAGE  2u:  PRINT 
1130  GOTO  1059 

1140  PRINT  :  PRINT  “LISTING  ENDED" 5  PRINT  CHRi  <12K  PR#  0 
1150  PRINT  “1  TO  LIST  AGAIN" 

1152  INPUT  “2  TO  RUN  PROGRAM  AGAIN" fG 

1153  IF  G  =  1  THEN  GOTO  1050 
1170  GOTO  8 

3000  PRINT  J  PRINT  M1  FOR  SINGLE  RECORD" 

3005  PRINT  "2  FOR  OVERLAYS  " 

3006  INPUT  "3  FOR  STAT  OPTIONS  “iG 

3007  IF  G  =  2  THEN  GOTO  6000 

3008  IF  G  =  3  THEN  GOTO  3500 

3012  INPUT  "FILE  NUMBER  ( 0-59  )J  "iZ 

3013  PRINT  5  PRINT  "ACCESSING  FILE  B">Z:  GOSUB  9500 

3040  IF  ER  =  1  THEN  PRINT  "FILE  EMPTY" :ER  =  05  GOTO  3012 
3090  PRINT  {  PRINT  "DATA  RETRIEVED":  PRINT 

388  am  =i  m 

3 m  m  ?  km.  ?oto  1050 

3110  XM  =  o:  GOSUB  7000 
3129  G  =  0:  GOSUB  8000 

3200  INPUT  Ai.*  IF  Aft  =  "P"  THEN  GOSUB  9000 
3210  CALL  62450:  HOME  5  TEXT 

3220  INPUT  "PRESS  RETURN  FOR  CYCLE  FREQUENCY  "?A#:  HGR2 

3240  XM  =  0:  GOSUB  7000 

3250  FOR  I  =  0  TO  (60  *  MINS)  -  1 

3260  FR  =  PEEK  (DSTART  +  2  +  <3  »  I)) 

3280  HCQLOR=  7J  HPLOT  I  +  10* 128  TO  I  +  10>128  -  (6.4  ft  FR) 

3290  NEXT  IJ  INPUT  Ai:  IF  Ai  =  “P"  THEN  GOSUB  9000 
3310  HOME  :  TEXT  J  CALL  624505  GOTO  1150 
3500  REM  ftSTAT  ROUTINE* 

3505  HOME 

3510  PRINT  :  PRINT  "(1)  SHOW  TRIAL  AND  MEAN*  SD" 

3530  HGR2  :  HCOLOR=  75XM  =  0:  GOSUB  7000 
3540  GOSUB  8000 

3550  PRINT  Di!"BLQAB  BSTAT  *A25561" 

3555  FOR  I  =  0  TO  (60  ft  MINS)  -  1 

3540  N  =  PEEK  ( SSTART  t<3*  I)):  IFN=0  THEN  GOTO  3610 
3565  XSUM  =  PEEK  ( SSTART  +14(3*1)) 

3570  X2SUM  =  PEEK  (SSTART  +  2  +  (3  ft  I)) 

3530  SBEV  =  10  ft  <  SQR  ((N  ft  X2SUM )  -  (XSUM  t  2)))  /  N 
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I 

1 


3585 

3590 

3595 

3600 

3605 

3610 

3615 

3621 

3622 
3640 
6000 
6020 
6025 
6040 
6060 

6070 

6071 
6075 

6080 

6085 

6100 

6110 

6150 

6155 

6160 

6170 

6190 

6230 

6235 

6240 

6245 

6260 

6270 

6280 

6290 

7000 

7010 

7020 

7030 

7040 

7050 

7055 

7060 

7061 

7062 

7063 

7064 

7065 

7066 

7067 
7070 
8000 
8001 
8002 
8003 
3010 
8012 
3013 

8014 

8015 
8020 
8021 
8022 

8023 

8024 
8030 


OViil 

8050 

8095 

3099 

9000 

9005 

9010 

9012 

901$ 

9020 

9500 


51  =  XBAR  4  SDEV!  IF  SI  >  100  THEN  SI  =  100 

52  =  XBAR  -  SDEV!  IF  S2  <  0  THEN  S2  =  0 

HPLOT  I  4  10r 128  -  ( 1.28  $  SI) 

HPLOT  I  4  10.123  -  (  1.28  *  S2) 

NEXT  I 

INPUT  A$:X  =  FRE  <0>:  IF  A*  =  "P"  THEN  GOSUB  9000 
CALL  624501  HONE  *.  TEXT 

PRINT  •  (LAST  FILE  ADDED*  "5  PEEK  (25561)."  )* 
PRINT  “  (THIS  FILE  IS  #",*Z>")" 

GOTO  3000 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


PRINT  "VARIABLES!  1  *  OVERALL  X" 


PRINT 

INPUT 

INPUT 

INPUT 


2  =  INSTANTANEOUS  V 

3  =  CYCLE  FREQUENCY 
1  5TR 

!UR 


'»G 


1  *TRi“  THROUGH  ‘.UR!  GOTO  61 


ist  file: 

2ND  FILE! 

PRINT  :  PRINT  "1  FUR  ONLY  THESE  TWO!" 

INPUT  “2  FOR  ALL  FILES  IN  BETWEEN." i ZR 
IF  ZR  =  2  THEN  PRINT  J  PRINT  “  FILES 
00 

PRINT  !  PRINT  "ACCESSING  VARIABLE  ",G 
PRINT  "  FILES  " »TR»“  L  "iUR 

7  a  TP*  npC|ir>  Qc,(\Ci 

IF  Er‘=  1  THEN  PRINT  "FILE  EHPTY" JER  =  0!  GOTO  6040 
HGR2  :  CALL  62450 
XM  =  0!  GOSUB  7000:  GOSUB  8000 
IF  ZR  =  2  THEN  Z  =  Z  4  1J  GOSUB  9500!  GOTO  6190 
2  =  UR!  GOSUB  9500 

IF  ER*=  1  THEN  TEXT  :  PRINT  "FILE  EMPTY" !ER  =  OJ  GOTO  6040 
GOSUB  3000 

IF  Z  <  UR  THEN  GOTO  6160 

INPUT  A$!  IF  A$  =  “P"  THEN  GOSUB  9000 

CALL  62450!  HOME  !  TEXT 

PRINT  "1  FOR  MORE  OVERLAYS." 

INPUT  "2  TO  EXIT.  "»G 
IF  G  =  1  THEN  GOTO  6000 
GOTO  8 
HCOLOR=  7 

HPLOT  5.0  TO  5.131  TO  270.131 
FOR  I  =  0  TO  10 

J  *  3!  IF  I  *  0  OR  I  *  5  OR  I  «  10  THEN  J  =  0 
HPLOT  J.12.S  *  I  TO  5.12.8  %  I 
NEXT  I 

IF  XM  =  1  THEN  GOTO  7070 
FOR  I  =  0  TO  3 
FOR  J  =  0  TO  5!K  =2 
IF  J  =  0  THEN  K  =  7 
IF  J  =  3  THEN  K  =  4 
QX  =  (I  *  60)  4  <J  *  10)  4  10 
HPLOT  QX. 131  TO  QX.131  4  K 
NEXT  J 
NEXT  I 
RETURN 

FOR  I  =  0  TO  (60  *  MINS)  -  1 
IF  G  =  1  THEN  GOTO  8010 

IF  G  =  2  THEN  GOTO  8020 

IF  G  =  3  THEN  GOTO  8030 

IF  G  =  0  THEN  HCOLOR*  7 

SY  =  128  -  (1.28  *  PEEK  ( DSTART  414(3*  I))) 

IF  SY  <  0  THEN  SY  =  0 
HPLOT  I  4  10, SY 
IF  G  >  0  THEN  GOTO  8050 

IF  G  =  0  THEN  HCOLOR*  7 

PY  *  128  -  (1.28  *  PEEK  (DSTART  4(3*  I))) 

IF  PY  <  0  THEN  PY  =  0 
HPLOT  I  4  10, PY 
IF  G  >  0  THEN  GOTO  8050 

IF  G  =  0  THEN  HCOLOR*  7 

FR  =  PEEK  (DSTART  424(3*1)) 

IF  G  =  0  THEN  HPLOT  I  4  10.128  TO  I  4  10.128  -  (FR)!  GOTO  8050 
HPLOT  I  4  10.128  -  (6.4  *  FR) 

G  —  G 
NEXT  I 
RETURN 

REM  SUBRT  TO  LINK  PRINTER  PROGRAM 
TEXT  !  PRINT  "STORING  ..." 

PRINT  Bt»"BSAVE  HRES,A16384,L8191,D1" 

PRINT  5  PRINT  "LOADING  GRAPHICS  PROGRAM... "!  PRINT 
PRINT  D$1 "RUN  PARA  PRINT" 

RETURN 

REM  SUBRT  TO  RETRIEVE  A  FILE 


PAGE  38 


1 

I 

1 

1 

I 


9501 

9505 

9510 

9511 

9512 

9513 

9514 

9515 

9516 

9517 
9513 

9519 

9520 
9580 
9300 

9801 

9802 

9803 

9804 
9810 
9815 
9820 

9330 

9860 

9866 

9867 
9368 
9870 
9880 


REM  NEEDS  MZM  FROM  MAIN  PROGRAM 
IF  PEEK  ( L'START  +  Z  )  =  0  THEN  ER  =  1 1  GOTO  9580 
PRINT  DtJ-BLOAD  B" »Z»" .A" JDSTART  -  10 
SBNUM  =  PEEK  (24576) 

CBN  =  PEEK  (24577) 

TL  =  PEEK  (24573) 

DAY  =  PEEK  (24579) 

MTH  =  PEEK  (24530) 

YEAR  =  PEEK  (24581) 

HOUR  =  PEEK  (24532) 

MNTS  =  PEEK  (24583) 

MINS  =  (  PEEK  (24534))  /  10 
MVOLTS  =  PEEK  (24585) 

RETURN 

REM  ID  DISPLAY  SUBROUTINE 

PRINT  {  PRINT  "I  TO  INSPECT":  PRINT  {  INPUT  "P  TO  PRINT:  "?A$ 

CP  =  C:  IF  A*  =  “P"  THEN  CP  =  1 
IF  CP  =  0  THEN  GOTO  9810 
PR*  i:  PRINT  CHR$  (9);"80Nu:  PRINT  CHR$  (30):  PRINT  CHR$  (01) 
PRINT  "S*  ".SBNUM*"  /  CONDITION  “ ‘ 


PRINT 
PRINT 
IF  MNTS 
JYEAR:  GOTO  9860 


:hour.*“:o";mnts;"  on 


PRINT  "AT 
PRINT 

PRINT  “LENGTH:  ".MINS."  MINUTES" 5  SPC(  3  )i "LEVEL? 
LTS" 

PRINT  " - " 

IF  CP  =  1  THEN  PRINT  :  PRINT  CHR$  (02)J  GOTO  9830 
INPUT  A$ 

RETURN 


CBNf "  /  TRIAL  " f TL 

ON  "JMTH»“/*»*DAY»V 

.‘mth}“/";day;"/"jyear 

.MVOLTS."  MICROVO 


I 

1 

I 

] 

l 

i 

l 

i 

l 

J 


a 
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f  TSYNTAX  ERROR 

.4  3PROGRAM  TO  LIST  SUMMARIZED  DATA 

_  YSYNTAX  ERROR 
f  3LIST 

1  1  D$  =  CHRi  (A) 

A  HIMEM*.  16384  JUST ART  =  13000 
■  5  DEF  FN  C(Q>  =  PEEK  (WLOC  +  0) 


6  MAR  =  5J  REM  MARGIN  SET  AT  5 
8  CALL  62450 J  TEXT  :  HOME  IX  = 


CALL  62450*.  TEXT  :  HOME  IX  =  FRE  <0)J  GOTO  3000 
))  /  100 

DO  PRINT  5  PRINT  "LISTINGS  OF  SUMMARIZED  DATA  BY  SUBTEXT t" 

D1  INPUT  “FIRST  SUBJECT  #.*"? 6.*  IF  (B  <  1)  OR  <B  >  27)  THEN  GOTO  3001 

D2  INPUT  “FINAL  SUBJECT  *{"fC:  IF  <C  <  B)  OR  (C  >  27)  THEN  GOTO  3002 

03  GOTO  9900 

DO  REM  FILE  HEADER  PRINT  ROUTINE 


PRINT  Bi?“PR*l" 
(02) 

FOR  I  =  F  TO  C 


PRINT  CHRi  (9)5“80N“*.  PRINT  CHRi  (30)J  PRINT 


FOR  K  =  1  TO  3 
PRINT  CHRi  (11):  REM 


?USTART5“ *D2" 


PRINT  CHRi  (11):  REM  VERTICAL  TAB 
NEXT  K 

PRINT  SPC<  MAR)* “SUBJECT  #"?IJ  PRINT 
PRINT  SPC(  MAR  )* "FILE  *  CONDITION 


PRINT  SPC(  MAR  )* “FILE 
E  LEVEL  X  ALPHA 
PRINT 

FOR  J  =  0  TO  31 
WLOC  =  USTART  +  ( 10  *  J ) 
IF  FN  C<0)  =  0  THEN  Gl 
PRINT 

IF  FN  C<6>  >  9  THEN  0 
PRINT  SPC<  . 


FREQ.1 


GOTO  9936 


TRIAL*  LENGTH 


IF  FN  C<6>  >  9  THEN  GOTO  9934 

PRINT  SPC<  MAR  +  1)?J?  SPC(  10)5  FN  C(0)?  SPC(  8)*  FN  C(l)*  SPC(  8) 
*<  FN  C(  2 ) )  /  10?  SPC<  4)?  FN  C(  3  )*“/“  $  FN  C(4)?V'80“5  SPC(  3)?  FN  C 
»  FN  C(6)f  SPC(  5)?  FN  C(7)J  SPC(  3)?  FN  C(8)?  SPC<  7)?(  FN  C 
<9j)  /  100 
GOTO  9935 

PRINT  SPC<  MAR  +  1)?J?  SPC(  10)5  FN  C(0)?  SPC(  8)*  FN  C(l)?  SPC<  8) 
?(  FN  C( 2 ) )  /  10?  SPC(  4)?  FN  C(3)?“/“?  FN  C(4)?'/'80“?  SPC<  3)?  FN  C 
<5 )?*:“?  FN  C<  6)?  SPC(  5)?  FN  C(7)?  SPC(  8)5  FN  C(8)?  SPC<  7) 5<  FN  C( 
9))  /  100 
NEXT  J 

PRINT  CHRi  (12) 

NEXT  I 

PRINT  Di? "PR*0* 

PRINT  CHRi  (7)5  CHRi  (7) 

PRINT  “LISTING  COMPLETED" 

INPUT  "M  FOR  MORE,  0  TO  QUIT"5A$ 

IF  Ai  =  "MM  THEN  GOTO  8 
END 


1 
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’SYNTAX  FRRDR 

IPRINT  PROGRAM  TO  PFRFORM  RFGRFSSION  $  CORRFI  ATTON  ON  SUMMARY  DATA 
0 

’SYNTAX  FRROR 
11  1ST 

1  D*  =  CHR*  (4) 

4  htmfm:  16334  histart  =  nooo 

5  CALL  62450:  PR  TNT  "I  OADING  CHARACTFRS" 

6  PRINT  CHR*  <  4  } : * P.!  HAD  A/N*A8197":  pnKF  ?3?,0!  POKF  233*32 
8  CALL  62450:  TFXT  :  HOMF  :X  =  FRF  <0)5  GOTO  3000 

3000  PRINT  "RFGRFSSION  6  CORRELATION INPUT  “SUBJECT  #:MfG 

3001  GOTO  9900 

7000  HCOl.OR*  7:  REM  AXES  SUBRT 
7002  HPLOT  75,0  TO  75,131  TO  259,131 
7004  FOR  I  =  0  TO  10 

7006  J  *  3:  IF  I  =  0  OR  I  *  5  OR  I  =  10  THEN  J  =  0 
7008  HPLOT  J  +  70,12,3  *  I  TO  75,12.8  *  I 
7010  NEXT  I 
7012  FOR  I  =  0  TO  12 

7014  J  =  2t  IF  I  =  0  OR  I  =  5  OR  I  =  10  THEN  J  =  5 
7016  OX  =  (14,976  *  I )  +  78 

7018  HPLOT  GX, 131  TO  OX, 131  +  J 

7019  NEXT  I 

7020  SCALE*  l:  ROT=  0 

7022  OR AW  18  AT  60,1311  DRAW  16  AT  55,671  DRAW  23  AT  60,67 
7024  DRAW  19  AT  43,7:  DRAW  16  AT  55,71  DRAW  .18  AT  60,7 

7026  DRAW  13  AT  76,1461  DRAW  16  AT  149,146:  DRAW  23  AT  154,146 

7028  DRAW  19  AT  220,146:  DRAW  16  AT  227,1461  DRAW  18  AT  232,146 
7030  RETURN 

9600  REM  CONVERSION  SUBROUTINE 

9602  IF  PEEK  ( WLOC )  =  0  THEN  ER  =  1J  GOTO  9610 

9604  FRSUM  =  (  PEEK  (WLOC  +  9))  /  100 

9606  PER  =  (  PEEK  (WLOC  +  8))  /  100 

9610  RETURN 

9900  REM  CALC  ( REG/CORR  )  ROUTINE 

9901  INPUT  “LARGE  OR  SMALL  GRAPH?" ii% 

9902  HGR2  1  GOSUB  7000 

9903  XS  =  01XQ  =  01 YS  =  01 YQ  =  0!XY  =  OJN  *  01  SCALE*  11  ROT*  0:  IF  L*  =  * 
S"  THEN  SCALE*  2 

9904  PRINT  DtPBLOAD  S" »G»" »A"  5USTART»",D2" JZ  =  0 

9905  FOR  I  =  C  TO  31 

9906  UL0C  =  USTART  +  ( 10  *  I ) 

9907  GOSUB  9600 :  IF  ER  =  1  THEN  GOTO  9921 

9908  XS  =  XS  +  FRSUM1YS  =  YS  +  PER  f 

9909  XY  =  XY  4  ( FRSUM  *  PER) 

9911  XQ  *  XQ  4  (FRSUM  t  2  )1YQ  =  YQ  4  (PER  t  2) 

9913  N  =  N  4  1 

9915  YP  =  128  -  (1.28  *  (  100  *  PER)) 

9916  XP  =  78  4  (149.76  *  FRSUM) 

9917  DRAW  66  AT  XP,YP 

9920  NEXT  I 

9921  BN  =  (N  *  XY)  -  (XS  *  YSKBB  =  (N  I  XQ)  -  (XS  t  2  )JER  *  0 
9923  BF  *  BN  /  BD 

9925  AF  =  (  YS  /  N)  -  ( BF  *  (XS  /  N>) 

9927  RD  =  SQR  ((N  %  XQ)  -  (XS  t  2))  *  SQR  ( <  N  *  YQ )  -  ( YS  t  2 ) ) 

9929  RF  =  BN  /  RD1TF  =  ( RF  *  SQR  (N  -  2))  /  SQR  (1  -  RF  t  2) 

9930  XI  =  OJYl  =  AF 

9931  IF  AF  <  -  .15  THEN  XI  =  (  -  .15  -  AF)  /  BFJY1  =  -  .15 

9932  IF  AF  >  1.0  THEN  XI  *  (1  -  AF)  /  BF1Y1  =  1 

9933  X2  =  1.2JY2  =  (BF  *  1.2)  4  AF1PY  *  Y2 

9934  IF  PY  >  1.0  THEN  X2  =  ( 1  -  AF )  /  BF5Y2  =  1 

9935  IF  PY  <  -  .15  THEN  X2  *  (  -  .15  -  AF)  /  BFJY2  =  -  .15 

9936  HPLOT  (149.76  *  XI )  4  78,128  -  <1.28  *  <100  *  Y1 ) )  TO  (149.76  *  X2)  4 
78,128  -  (1.28  *  (100  *  Y2)> 

9937  PRINT  CHR*  (7),  CHR*  (7)J  INPUT  •  *1A* 

9939  TEXT  :  HOME  :  PRINT  "STORING  .  .  ." 

9940  PRINT  B*5"BSAVE  HRES,A16384,L8191,D1“ 5  PRINT  "STORED" 

9941  PRINT  B*i"PR*rJ  PRINT  CHR*  (  9  )i”S0N"  5  PRINT  CHR*  <  02  )t  PRINT  CHR* 
( 30 ) 

9942  PRINT  CHR*  (11):  PRINT  CHR*  (11):  PRINT  CHR*  (11)!  PRINT  CHR*  <1 
n:  PRINT  SPC<  22  )* "SUBJECT  *"}G!  PRINT 

9943  PRINT  5  PRINT  SPC<  24)f“X  =  FREQUENCY  OF  STATE  CHANGE  (FULL  CYCLE  I 
N  HZ )  i  PRINT 

9944  PRINT  3PC<  24)i"Y  =  PROPORTION  OF  TIME  IN  ALPHA".*  PRINT 

9946  IF  AF  <  0  GOTO  9948 

9947  PRINT  3PC(  26 ^"REGRESSION:  Y  *  “ J  BF  * " X  +  "JAFJ  GOTO  9949 
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l 


J 

I 

3 

3 

! 

3 

1 


9948 

9949 

9950 

VS 

9955 

9957 


PRINT 
PRINT 
PRINT 
40  )#  *N 
PRINT 
PRINT 
PRINT 
END 


SPC<  _ 
t  PRINT 
:  PRINT 
=  "  »N 
PRINT 


26  )* "REGRESSION*  Y  =  “»BF; 
SPCI  26)” CORRELATION:  - 
SPC(  - 


40  )i 


=  '* » TF I" 


: LOADING ' GRAPH IcS$PRGGRAN. 
0$ »  RUN  PARA  PRINT  2 fDl" 


■;af 

R  =  " 5  RF 
DF  =  "fN 


2t  PRINT  5  PRINT 


SPCI 


3PR40 

3 


PACE  U2 


’SYNTAX  ERROR 

3PRINT  1-WAY  ANOVAn  WITHIN  SUBJECT 
1  0 

3LIST 

1  Dt  =  CHR*  (4) 

3  BEE  FN  ASN(  X >  =  ATN  <X  /  SQR  (  -  X  *  X  +  1 ) ) 

5  BIN  SS(4)J  DIN  MSI  4)1  DIN  ME(2>!  DIM  DF2<4)!  DIM  F(  2 ) 

10  DIM  MAD(4>16>!  DIM  BUF<4»16) 

12  DIM  ST(16>!  DIM  CT(4> 

14  HONE 

15  PRINT  "REPEATED  MEASURES  ANQVA  *  * 

17  PRINT  “ONE-WAY *  WITH IN-SUBJECT" 

20  PRINT  1  INPUT  “4  OF  LEVELS  OF  VARIABLE!  "JV1 
22  INPUT  “VARIABLE  NAME!  "JAN* 

28  IF  (VI  >  4)  THEN  PRINT  "TOO  MANY  VARIABLE  LEVELS!  MUST  BE  <5" ! 


29 

30 

31 

32 

35 

36 
38 
40 
42 
46 
48 
50 


20 


■»G$ 


OJAS  =  OJTT  =  05XQ  =  0 


INPUT  ”*  OF  SUBJECTS!  "JSN 
INPUT  "DEPENDENT  VARIABLE  NAMF5 
PRINT  !  PRINT  "INPUT  DATA" 

FOR  I  =  1  TO  VI 
HOME 

FOR  K  =  1  TO  SN 
PRINT  "A" ? I J" >  SJJK 
INPUT  “VALUE  =  “ JBUF(IfK) 

NEXT  K 

NEXT  I!  HOME 

PRINT  "DATA  MATRIX  COMPLETE"!  GOSUB  400 
PRINT  "WORKING  ON  ANALYSIS" 

FOR  I  =  1  TO  4 

54  SS(I)  =  0!MS(  I )  =  0!CT(I)  =  0 

56  NEXT  I 

57  AO  =  OJTQ  =  OJSQ  =  0!SA  =  0!SB 

58  TS  =  0 

60  FOR  I  =  1  TO  16 

61  ST<  I )  =  0 

62  NEXT  I 

67  FOR  I  =  1  TO  VI 

69  FOR  K  =  1  TO  SN 

70  TT  =  TT  +  MAB< I»K ) 

71  CT(  I )  =  CT(  I )  +  MAD(  I»K  )!ST(  K )  =  ST(K)  t  MAD(I»K) 
74  NEXT  K!  NEXT  I 

80  DF%( 1 )  =  VI  -  1 ! BFX( 2 )  =  SN  -  1 

81  DF%(  3 )  =  DF%(  1 )  *  BFX(2)!BFX(4)  =  <  VI  *  SN  >  -  1 

82  T9  =  TT  t  2 

84  FOR  I  =  1  TO  VI 

85  FOR  K  =  1  TO  SN 

86  XQ  =  XQ  +  MAD(IfK)  t  2 

88  NEXT  K!  NEXT  I 

89  FOR  I  *  1  TO  VI 

90  AG  =  AQ  +  CT< I )  t  2 

91  NEXT  I 

92  FOR  I  =  1  TO  SN 

93  SO  =  SQ  +  ST(I )  f  2 

94  NEXT  I 

96  SA  =  (AQ  /  SN) 

98  SB  =  (SQ  /  VI) 

100  AS  =  XQ 

102  TS  =  (TT  t  2)  /  (VI  *  SN) 

104  SS( 1 )  =  SA  -  TS 

106  SS(2)  =  SB  -  TS 

108  SS( 3 )  =  AS  -  SA  -  SB  +  TS 

110  SS<4)  =  AS  -  TS 

135  ..EQR.I  ■  1  TO  3 
* 


1 


I 


5(1)  /  DFX(  I ) 


134  NEXT 

140  F(  1 )  s  MS(  1 )  /  MS(  3 )!F( 2 )  =  MS(2)  /  MS(3) 

144  FX  =  100  *  F(  1  )!F(  1 )  =  FX  /  100!FX  =  100  *  F(2)5F(2)  =  FX  /  100 
146  HONE  !  PRINT  "ANOVA  SUMMARY  TABLE!  "JAN*J"  BY  SUBJECTS" 

PRrNT  11 _ — _ — _ _ _ « 

PRINT  "SOURCE"  J  SPC<  2)»"S.S."»  SPC(  6)J,,D.F."J  SPC<  4)J"M.S.HJ 
7  )i  ”F” 

PRINT  " - » 

154  DATA  *  A  S'  S  AXS  1 

156  RESTORE 

160  FOR  I  =  1  TO  3 


143 

150 


GOTO 


SPC< 


162  READ  OIJ  IF  (<T  =  j)  hr  <r  s  o))  THEN  GOTO  166 

164  PRINT  5  PRINT  3$J  SRC?  l)»3S(t)»  TABC  19  >»DFZC  I  >*  TAB(  23  )JMS<  I  )J  GOTO 

168 

166  PRINT  |  PRINT  Gif  SPC(  l)fSS(I)f  TABC  19)fDFX<I)f  TABC  23)fMS(I)»  TABC 
3j  ;>F(  I ) 

168  NEXT  I 

169  PRINT  " - ’ 

170  PRINT  "TOTAL"  f  SPCC  lHSSUW  SPCC  2)*DF%C4) 

172  INPUT  ZU  RESTORE 

171  TiJptif  “HiST!  :'Pi>v?  t  V / it  1  »*7* 

176  IF  Z4  =  * Y*  THEN  PRINT  D«f"PR*l":  PRINT  C.  (9)fM80N“:  PRINT  CHRI 

(30)1  PRINT  CHRI  (02)5  GOTO  180 

177  HONE  l  INPUT  "C  TO  CONTINUE*  N  FOR  NEW  DATA  At 

178  IF  At  =  “C*  THEN  GOTO  48 

179  HONE  t  GOTO  15 

180  PRINT  "ANGVA  SUGARY  TABLE  I  "5ANt*“  BY  SUBJECTS' 

182  PRINT  " - " 

184  PRINT  "SOURCE"}  SPC(  2)*MS.S.M}  SPC(  6)fHD.F."f  SPCC  4)f“H.S."f  SPCC 
7  )*"F" 

186  PRINT  ■ - " 

188  FCR  I  =  1  TO  3 

190  READ  at:  IF  CCI  *  1)  OR  Cl  *  2>)  THEN  GOTO  194 

192  PRINT  t  PRINT  Qtf  SPCC  1)*SS(I)»  SPCC  3 >?DF%<  I  >*  SPC(  2 )*«S(  I  )J  GOTO 

195 

194  PRINT  1  PRINT  Qtf  SPCC  l)fSS(I)f  SPCC  3)fDFZ(I)f  SPC(  2)fMS(I)f  SPC( 

2  )»F(  I ) 

195  NEXT  I 

196  PRINT  " - " 

19?  EgtMI  "Ictal";  spcc  ittssu);  spcc  jlurfzu):  gosub  300 

198  PRINT  CHRI  (12)1  PRINT  Btf'PRlO  *  GOtO  1/2 
300  REN  PRINT  NEANS 

302  PRINT  PRINT  *.  PRINT  "CELL  MEANS!":  PRINT 
304  FOR  I  =  1  TO  Vi 

308  PRINT  "A"* If"  =  ";CT(I>  /  SN 

309  PRINT 

310  NEXT  I 

314  IF*TFRM  >  4  THEN  A$  =  "l/“ 

315  IF  TFRN  =  4  THEN  A$  =  "LOG" 

316  IF  TFRN  =  3  THEN  At  =  "ARCSINE" 

317  IF  TFRN  =  2  THEN  A$  =  "SQUARE  ROOT" 

318  PRINT  "DEPENDENT  VARIABLE  =  ‘SGI 

319  IF  TFRN  >  1  THEN  PRINT  "TRANSFORMED  BY  " ;A$;"(X+“ fCNf" )" 

3?q  RETURN 

400  REN  PICK  TRANSFORMATION 

401  TFRN  =  0:  PRINT  “TRANSFORMATION:":  PRINT 

402  PRINT  "1  =  NONE" 

404  PRINT  ”2  =  SQUARE  ROOT" 

406  PRINT  "3  =  ARCSINE" 

408  PRINT  “4  =  L0G<  X )* 

409  INPUT  "5  =  l/X"fTFRM 

410  CN  =  0 

412  PRINT  :  INPUT  "ADD/SUBTRACT  CONSTANT:  "fCN 

413  PRINT  "TRANSFORMING" 

414  FOR  I  =  1  TQ  VI 

418  FOR  K  =  1  TQ  SN 

419  U  =  BUF(IfK) 

420  IF  TFRM  >  4  THEN  T  =  1  /  <  U  +  CN ) 

422  IF  TFRM  =  4  THEN  T  =  LOG  (U  +  CN) 

424  IF  TFRM  =  3  THEN  T  =  FN  ASN(  W  +  CN) 

426  IF  TFRM  =  2  THEN  T  =  SQR  (U  +  CN  ) 

428  IF  TFRM  <  2  THEN  T  =  U 

430  KADCI*fO  =  T 
432  NEXT  K:  NEXT  I 
434  RETURN 
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]PR*0 


1 SYNTAX  ERROR 

'PRINT  2-UAY  AMOVA*  UITHIN-SUBJECT 
2  0 

31  1ST 

1  ?;  -  CHRI  (4) 

3  PEF  FN  ASN<  X )  =  ATN  <  X  /  SQR  (  -  X  *  X  +  J  > ) 

5  PIM  SS(8):  DIM  MS(S)J  DIM  MEU.4K  DIM  DFX(8> 

10  DIM  KAD<4*4*16>:  DIM  SA(4*16).*  DIM  SB(4*16)J  DIM  BUFU.4.16) 

14  HOME 

15  PRINT  "REPEATED  MEASURES  ANOVAJ" 

17  PRINT  "TUO  VARIABLES*  UITHIN-SUBJECT" 

20  PRINT  1  INPUT  "*  OF  LEVELS*  VARIABLE  15 


22  INPUT  “VARIABLE  NAME;  SAN* 


24 

26 

28 

29 

30 

31 

32 

34 

35 

M 

39 

40 

41 

42 
44 
46 
48 
50 
52 


65 

67 

68 
69 


PRINT 


*V1 
1  *V2 


:'5»  : 

INPUT 

INPUT 

PRINT 


PRINT  "TOO  MANY  VARIABLE  LEVELS?  MUST  BE 


*Gt 


INPUT  "*  OF  LEVELS*  VARIABLE  25 
INPUT  "VARIABLE  NAME?  "*BN$ 

IF  (VI  >  4)  OR  <V2  >  4)  THEN 
GOTO  20 

*4  OF  SUBJECTS:  SSN 
'DEPENDENT  VARIABLE.* 

5  PRINT  "INPUT  DATA" 

FOR  I  *  1  TO  VI 
FOR  J  =  1  TO  V2 
HOME 

FOR  K  =  1  TO  SN 
PRINT  "A" *1 5" *  B" * J*  *  S  *K 
REM  GOTO  41  SOMETIMES 
INPUT  "VALUE  =  " *BUF(  I *J*K ) 

REM  READ  MAD(I*J*K)  SOMETIMES... 

NEXT  K 
NEXT  J 

NEXT  1 5  HOME 

PRINT  "DATA  MATRIX  COMPLETE":  GOSUB  400 
PRINT  " UORKINS  ON  ANALYSIS" 

FOR  I  =  1  TO  8 
54  SS(I  >  =  05MS(  I ) 

56  NEXT  I 

57  FOR  I  =  1  TO  4 

58  FOR  J  =  1  TO  12 

59  SA(  I  *  J )  =  0:SB(I*J)  =  0 

60  NEXT  j:  NEXT  I 

61  FOR  I  =  1  TO  4 

62  FDR  J  =  1  TO  4 

63  ME(  I*J  )  =  0 

64  NEXT  J{  NEXT  I 

■  t  =  oja  =  o:s  =  o:as  =  o:b  =  o:bs  =  o:ab  =  o:axbxs  =  o 

FOR  I  =  1  TO  VI 
FOR  J  =  1  TO  V2 
..  FOR  K  =  1  TO  SN 

70  T  =  T  +  MAD(I*J.K):AXBXS  =  AXBXS  +  ( MAD(  I*J*K ) )  t  2 

71  SA(  I  *  K )  =  SA(I.K)  +  MAD(  I*  J*K ) 

72  SB(J*K)  =  SB(  J*K )  +  MAD(I*J*K) 

73  ME(I.J)  =  ME(I»J)  +  MAD(  I*J*K ) 

74  NEXT  KJ  NEXT  Jl  NEXT  I 
REM  GOSUB  200  SOMETIMES... 


=  0 


T  =  (T  t  2)  /  (VI  *  V2  *  SN) 

80  FOR  I  =  1  TO  VI 

81  K  =  0 

82  FOR  J  =  1  TO  V2 

84  K  =  K  +  ME(  I*J) 

85  AB  *  AB  +  (ME(IfJ))  t 

86  NEXT  J  J A  =  A  +  K  t  2 

88  NEXT  ISA  =  A  /  (V2  * 

90  K  =  0 

92  FOR  J  *  1  TO  V2 

93  K  =  0 

94  FOR  I  =  1  TO  VI 
96  K  =  K  +  ME(  I  *  J ) 

98  NEXT  IJB  *  B  +  K  t  2 

99  NEXT  JJB  =  B  /  (VI  *  SN) 

FOR  I  =  1  TO  SN 


SNKAB  =  AB  /  SN 


101  K  =  0 


AS  + 


FOR  J  =  1  TO  VI 
K  »  K  +  SA(  J  *  I > : As 
NEXT  JJS  =  S  +  K  t  25K  * 
FOR  J  *  1  TO  V2 


(  SA(  J » I ) ) 
0 


t  2 


■  no 

1  112 
*  4  «  i 


„  ffl 

1  234 

*  mi 


BS  =  BS  4  (SB(  J.I ))  t  2 
NEXT  J 

NEXT  I  IAS  =  AS  /  V2?BS  =  BS  /  VI  IS  =  S  /  ( VI  *  V?  ) 

SS(1)  =  A  -  T?SS<2)  =  S  -  T?SS<3)  =  as-a-s  +  t 
SS(4)  =  B  -  T?SS<5)  =  BS  -  B  -  S  4  T?SS<6)  =  AB-A-B4T 
SS( 7 )  =  AXBXS  -AB-AS-BS4A4B4S  -  T?SS(8>  =  AXBXS  -  T 
DFX(l)  =  VI  -  1?»FX(2>  =  SN  -  1?DFX<3)  =  (VI  -  1)  *  (SN  -  1) 

BFX( 4 )  =  V2  -  1?DFX<5>  =  <V2  -  1)  *  (SN  -  1 ) J DFZ<  6 )  =  (VI  -  1)  *  ( V2  - 

DFX(  7 )  =  (VI  -  1)  *  (V2  -  1)  *  (SN  -  1>?HFX(8>  =  ( VI  *  V2  *  SN >  -  1 
FOR  I  =  1  TO  7 
MS<  I  >  =  SS(  I )  /  PFX(  I ) 

NEXT  I 

rl  =  MS(  1 )  /  MS(3)?F2  =  HS(4)  /  MS(5) 

F3  =  MS(  6 )  /  MS(  7  ) 

FX  =  100  *  F1JF1  =  FX  /  100JFX  =  100  *  F2JF2  =  FX  /  1005FX  =  100  *  F3J 
F3  =  FX  /  100 

F(  1 )  =  F1?F(4)  =  F2?F(6)  =  F3 
HOME  J  PRINT  "ANOVA  SUMMARY  TABLE?  " »AN*»"  BY  "JBN* 

PRINT  * - ■ 

PRINT  "SOURCE".  SPC(  2 ) » “  S .  S .  5  SPC(  6>»"D.F."»  SPC(  4  )f  "M.S. "  S  SPC( 

7  )>"F" 

PRINT  “ - ■ 

DATA  "  A  "»"  S  "»"  AXS  "»"  B  l»"  BXS  “»"  AXB  ".“AXBXS" 

RESTORE 

FOR  I  =  1  TO  7 

READ  3*?  IF  ((I  =  1)  OR  (I  =  4)  OR  (I  =  6))  THEN  GOTO  166 
PRINT  ?  PRINT  3*5  SPC(  l)fSS(I)?  TAB<  19  >$OFZ(  I  )5  TAB(  23)»HS(I)?  GOTO 
168 

PRINT  ?  PRINT  3*5  SPC(  l)iSS(I)}  TAB(  19)»DFX(I>.*  TAB(  23)»MS(I)5  TAB( 
35  )»F(  I ) 

NEXT  I 

PRINT  " - " 

PRINT  “TOTAL"?  SPC(  1)»SS(8>;  SPC(  2)fDFX(8) 

INPUT  Z*?  RESTORE 
INPUT  °HARti  COPY?  ( Y/N )  ” 

if  z*  =  “Y*  then'  p'rint  d$;,pr*i“?  print  chr*  (9);b80N“?  print  chr* 

(30)?  PRINT  CHRi  (02)?  GOTO  180 

HOME  ?  INPUT  "C  TO  CONTINUE.  N  FOR  NEW  BATA"5A* 

IF  A*  =  "C“  THEN  GOTO  48 
HOME  ?  GOTO  15 

PRINT  "ANOVA  SUMMARY  TABLE?  "JAN*."  BY  " »BN* 

PRINT  " - "  • 

PRINT  “SOURCE".*  SPC(  2)»"S,S.";  SPC(  6)*"D.F. “  J  SPC(  4);tM.S.“;  SPC( 
7WF" 

PRINT  " - " 

FOR  I  =  1  TO  7 

READ  3*?  IF  ((1  =  1)  OR  (1  =  4)  OR  (I  =  6))  THEN  GOTO  194 

PRINT  ?  PRINT  3*J  SPC(  DJSSd).  SPC(  3)»DFX(I)»  SPC(  2).‘MS(I)?  GOTO 

195 

PRINT  ?  PRINT  3*;  SPC(  1).*SS(I).*  SPC(  3)fDFX(I)»  SPC(  2).‘MS(I)f  SPC( 

2)»F(I) 

NEXT  I 

PRINT  " - " 

PRINT  "TOTAL".*  SPC(  1>.SS(8)5  SPC(  2)»*DFX(8)?  GOSUB  300 

PRINT  CHR*  (12)?  PRINT  D*»"PR*0"?  GOTO  172 

HOME 

PRINT  "ABS  MATRIX" 

FOR  I  =  1  TO  VI 
FOR  J  =  1  TO  V2 
FOR  K  =  1  TO  SN 
PRINT  K»  SPC(  4)»MAD(  I.J.K) 

NEXT  K?  INPUT  Z*?  NEXT  J?  NEXT  I 
INPUT  Z$ 

HOME  ?  PRINT  "AS  MATRIX" 

FOR  I  =  1  TO  SN 

PRINT  I,*  SPC(  4 )»SA(  1 .1 )»  SPC(  4).‘SA(2.I) 

NEXT  I?  INPUT  Z* 

HOME  ?  PRINT  "BS  MATRIX" 

FOR  I  =  1  TO  SN 

PRINT  I.‘  SPC(  3)»SB(1.I)»  SPC(  3)JSB(2.I)?  SPC(  3 )»SB(  3.1 )»  SPC(  3)»S 
B(4.I ) 

NEXT  I?  INPUT  Z* 

HOME  ?  PRINT  “AB  MATRIX" 

FOR  1=1  TO  VI 

PRINT  I;  SPC(  3)JME(I»1);  SPC(  3)5ME(  1.2  )J  SPC(  3)5ME(I.3).  SPC(  3)JM 

N^XT4^?  INPUT  Z* 

RETURN 

REM  PRINT  MEANS 
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1 

1 

1 

1 

1 

1 

* 

\ 


302 

304 

306 

308 

309 

310 

312 

314 

315 

316 

317 

313 

319 

320 

400 

401 

402 
404 
406 

408 

409 

410 

412 

413 

414 
416 

418 

419 

420 
422 
424 
426 
428- 
430 
432 
434 


PRINT  {  PRINT  ?  PRINT  'CELL  MEANS? *  5  PRINT 
FOR  1=1  TO  VI 
FOR  J  =  1  TO  V2 

PRINT  "A“»I»“  =  " J HE( I  * J )  /  SN 

PRINT 

NEXT  J?  NEXT  I 
PRINT 

IF  TFRM  >  4  THEN  A*  =  "1/" 

IF  TFRM  =  4  THEN  A$  =  "LOG" 

IF  TFRM  =  3  THEN  AS  =  ■ARCSINE" 

IF  TFRM  =  2  THEN  A$  =  ’SQUARE  ROOT" 

PRINT  "DEPENDENT  VARIABLE  =  "JG$ 

RETURN**  >  1  THEN  PRINT  “TRANSFORMED  BY  " iA$»"<X4" JCNi"  )" 

REM  PICK  TRANSFORMATION 
TFRM  =  0:  PRINT  "TRANSFORMATION:*:  PRINT 

1  =  NONE" 

2  =  SQUARE  ROOT" 

3  =  ARCSINE" 
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PRINT 
PRINT 
PRINT 
PRINT  "4  =  LOG<  X  )" 


INPUT 
CN  =  0 
PRINT 
PRINT 
FOR  I 
FOR  J 


1/X” i TFRM 

:  INPUT  "ADD/SUBTRACT  CONSTANT?  "»CN 
'TRANSFORMING" 


1  TO  VI 
=  1  TO  V2 
FOR  K  =  1  TO  SN 
U  =  BUF(IfJ»K) 


IF  TFRM  >4  THEN 
IF  TFRM  =  4  THEN 
IF  TFRM  =  3  THEN 
IF  TFRM  =  2 
IF  TFRM  <  2 
HAD(IfJtK)  = 

NEXT  K?  NEXT  J?  NEXT  I 
RETURN 


THEN  T  = 
THEN  T  =  U 
T 


1  /  (U  +  CN) 

LOG  (U  +  CN) 

FN  ASN(  U  +  CN) 
SQR  (V  +  CN) 


! 


ANOVA  S.'v„iASV  table:  training  by  subjects 
SOURCE  S.S.  D.F.  H*S.  ‘  F 


A  561 . 124992  1  561.124992  3.72 

S  11223  15  748.533332  4.96 

AXS  2260.87501  15  150.725001 

TOTAL  14050  31 

CELL  means: 

A1  =  23.8125 
A2  =  32.1875 

DEPENDENT  VARIABLE  =  %  ALPHA 


* 


Table  2.  Summary  of  analysis  of  variance  of  log  cycle  frequency  categorical 
by  training  exposure. 


ANOVA 

SUMMARY  TABLEt 

TRAINING  BY 

SUBJECTS 

SOURCE 

:  s.s. 

P.F 

.  K.S. 

F 

A 

.193226069 

1 

.193226069 

1.37 

S 

34.329386 

15 

2.28862573 

16.23 

AXS 

2.11394075 

15 

,140929383 

TOTAL 

36.6365528 

31 

CELL  means: 

A1  =  -.956967284 
A2  =  -1.11238046 

DEPENDENT  VARIABLE  =  CYCLE  FREQUENCY 
TRANSFORMED  BY  L0G(  X+0  ) 
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Table  3 


Summary  of  analysis  of  variance  of  %  alpha  categorical  by  feedback 
conditions  and  eye  condition. 


ANOVA 

r  4 

r  r  r  u  t*.A.  ...  n  » 

r  *  r  S 

snisar. 

r  S  *  3  * 

*  p  4  n  * 

r 

w 

S+U,  09799  1 

ivi<4  4  VT  T  9  7  1 

s 

5773.09993 

? 

•44  j 

AXS 

702.400009 

r* 

7 

73,0444^ 

B 

3534.39999 

5 

37.34  ♦  39999 

T  «  I 
r  «  »  *.* 

BXS 

1+72.10001 

? 

j  Aif » 

AXB 

448.900009 

1 

448,900009 

8*33 

AXBXS 

484.599991 

9 

B444434 

TOTAL  12919.6  39 


CELL  MEANS! 
A1  B1  =  +5.7 
A1  B2  =  20 .2 
A2  B1  =  31.9 
A2  B2  =  19.3 


DEPENDENT  VARIABLE  =  RAW  %  $  FnA 
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Table  4 


Summary  of  analysis  of  variance  of  baseline-scaled  %  alpha  cate¬ 
gorical  by  feedback  condition  and  eye  condition. 


ANOVA 

SUMMARY  TABLE.* 

FFFBhCK  by 

FYFS 

SOURCE  S.S. 

D.E 

.  H.S. 

F 

A 

577.600001 

1 

577.600001 

7.28 

S 

3528.1  9 

392 

.011111 

AXS 

713.900003 

9 

79.3222226 

B 

32.4000001 

1 

32.4000001 

.17 

BXS 

1621.1  9 

180 

,122222 

AXB 

384.399999 

1 

334.399999 

7.79 

AXBXS 

444.100001 

9 

49.3444445 

TOTAL  7301 .6  3? 

CELL  means: 

At  B1  =  9.4 
At  B2  =  5 
A2  B1  =  -4.4 
A?  B2  =  3.6 

BFPFNPENT  VARIABLE  =  BASELINE-SCALED  X  ALPHA 
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Table  5 


Summary  of  analysis  of  variance  of  log  baseline-scaled  %  alpha 
categorical  by  feedback  condition  and  eye  condition. 


SUMMARY  TAB* F! 

FFEBACK  BY 

F  tr  -i 

SOURCE  s.s* 

P.F 

.  M.S. 

F 

A 

1.13943768 

1 

1.17943768 

4.73 

s 

6.97473693 

9 

.77497077 

AXS 

2.16672468 

9 

.760747187 

B 

.54594624 

1 

.54594624  1 

:  .04 

BX5 

4.77.13093 

9 

,574589973 

AXB 

1.16587758 

1 

1.16587758 

5.03 

AXBXS 

2.08203995 

9 

.73133777? 

TOTAL 

18.7960725 

39 

CELL  MEANS? 

A1  B1  =  3  •  5907?-*  i  3 
A1  B2  =  3*48292954 
A2  B1  =  2.91171905 
A2  B2  =  3*4868235 


DEPENDENT  VARIABLE  =  BASELINE-SCALED  2  Al  PHA 
TRANSFORMED  BY  l.0G(X+30) 
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Table  6 


Summary  of  analysis  of  variance  of  reciprocal  %  alpha  by  feedback 
mode  and  feedback  polarity. 


riivuUn 

Tm rvl  Fj 

i'.nriF  BY  POLARITY 

SOURCE  S.S.  D.F 

• 

r:,s,  F 

A 

1.2434077.7F-03 

1 

) .24340737E-03 

7.15 

S 

1.99496974F-03 

5 

3.98993948F-04 

AXS 

2.88833869F-03 

5 

5.7766773RF-04 

B 

8.08664699F-05 

1 

3.08664699E-05 

.93 

BXS 

4.33332789E-04 

5 

•5,66665578F-05 

AXB 

2.73883998F-05 

1 

2.73883998F-05 

.41 

AXBXS 

3.32092073E-04 

5 

6.64134147E-05 

TOTAL 

7.00039553F-03 

23 

CELL  MFANS*. 

A1  B1  =  ,0223002561 
A1  B2  =  ,028)079846 
A2  B1  =  ,0388324235 
A2  B2  =  .040367105 


PEPENDFNT  VARIABLE  =  RAW  %  ALPHA 
TRANSFORMED  BY  l/CX+0) 
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Table  7 


Summary  of  analysis  of  variance  of  arcsine  cycle  frequency  by  feed¬ 
back  mode  and  feedback  polarity. 


vr* 

p?;;:r:ASV  tab: 

f: 

KOBE  BY  POLARITY 

soiiRcr 

7  S.S. 

D.F. 

,  H .  S .  F 

A 

.0863848105 

1 

.0883848105  .74 

s 

.961924925 

5 

.192334985 

AXS 

.579917073 

5 

.115983415 

B 

.0231894975 

1 

.0231894925  3.26 

BXS 

.0355340988 

5 

7.10681975F-03 

AXB 

.0278338417 

1 

.0278338417  7.74 

AXBXS 

.0179584473 

5 

3.59168947E-03 

TOTAL  1.73274269  23 

CELL  means: 

A1  B1  =  .977930006 
A1  F2  =  .847651548 
A2  B1  =  .789830559 
A2  B2  =  .795772054 

DEPENDENT  VARIABLE  =  RAW  CYCLE  FREQUENCY 
TRANSFORMED  BY  ARCSINE! X40  ) 


Summary  of  analysis  of  variance  of  square-root  baseline-scaled  t 
alpha  by  feedback  mode  and  feedback  polarity. 


ANOVA 

SUMMARY  TADL 

e: 

MODE  BY  POLARITY 

SOURCE  S.S. 

P.F 

.  M.S* 

F 

A 

1.47294509 

1 

1.47294509 

.11 

S 

16.1641321 

3 

5.38804404 

AXS 

39.1536011 

3 

13.0512003 

B 

3.9291S3S4 

1 

3.92918384 

1.01 

BXS 

11.6306747 

3 

3.87689157 

AXB 

1.30974019 

1 

1.30974019 

3.3 

AXBXS 

1.18736744 

3 

.395789146 

TOTAL 

74.8476446 

15 

CELL  MEANS* 

A1  B1  =  5*95624375 
Ai  B2  =  4.39291587 
A2  B1  =  4.77719956 
A2  B2  =  4.35831059 

DEPENDENT  VARIABLE  =  BASELINE-SCALED  X  ALPHA 
TRANSFORMED  BY  SQUARE  R00T(X+25) 
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Table  9 


Summary  of  analysis  of  variance  of  arcsine  baseline-scaled  cycle 
frequency  by  feedback  mode  and  feedback  polarity. 


ANOVA 

SUMMARY  TABLE? 

MODE  BY  POLARITY 

SOURCE  S.S.  D.F 

.  H.S.  F 

A 

•038361 1631  1 

.0388611681  .38 

S 

.0183447062  3 

6.28156875E-03 

AXS 

.30188102  3 

.100627007 

B 

1. 962251 26E-03 

1  1.96225126E-03 

BXS 

4.09634924E-03 

3  1.36544975E-03 

AXB 

.0118594922  1 

.0118594922  5.1 

AXBXS 

6.9634621E-03 

3  2.321 15404E-03 

TOTAL 

.38446845  15 

cell  means: 

A1  B1 

=  .323492707 

A1  B2 

=  .246893408 

A2  B1 

=  .170475379 

A2  B2 

=  .202777873 

DEPENDENT  VARIABLE  =  BASELINE-SCALED  CYCLE  FREQUENCY 
TRANSFORMED  BY  ARCSINE< X+.29 ) 


APPENDIX  D: 

CORRELATIONS  AND  REGRESSION 
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ALL  SUBJECTS 

X  =  FREQUENCY  OF  STATE  CHANCE  (FULL  CYCLE  IN  HZ) 
Y  =  PROPORTION  OF  TIME  IN  ALPHA 

REGRESSION:  Y  =  .549149122X  -.0115604562 

correlation:  r  =  .710700073 

T  =  20.3562678  DF  =  406 
N  =  408 


Figure  12. 


Correlation  and  regression:  scatter  plot  of  %  alpha  by  cycle  fre¬ 
quency  across  all  trials  and  subjects. 
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Alt  SUBJFCTSf  CONDITIONS  1»  3»  5t  7.  8 

X  =  FREQUENCY  OF  STATE  CHANGE  (FULL  CYCLE  IN  HZ) 
Y  =  PROPORTION  OF  TINE  IN  ALPHA 

REGRESSION:  Y  =  .503360182X  +  .0583454899 

CORRELATION:  R  =  .59&37849A 

T  =  11, 6792026  DF  =  ?45 

N  =  ?4  7 


0  .5  1.0 


Figure  13.  Correlation  and  regression:  scatter  plot  of  %  alpha  by  cycle  fre¬ 
quency  across  all  subjects  for  eyes-closed  rials. 


AH  StlB-TPHTS,  CONDITIONS  ?,  4,  6»  9.  10 


X  =  FRFOUFNCY  OF  STATE  CHANGE  (FULL  CYCLE  IN  H7 ) 
Y  =  PROPORTION  OF  TIMF  IN  Al  PHA 

REGRESSION:  Y  =  .436940376X  -.0336601179 

CORRELATION:  R  =  .S77674437 

T  =  23.0920964  OF  =  159 


N  =  161 


Figure  14.  Correlation  and  regression:  scatter  plot  of  t  alpha  by  cycle  fre¬ 
quency  across  all  subjects  for  eyes-open  trials. 


